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We describe in detail how the different components of a multi-phase transport (AMPT) model, 
that uses the Heavy Ion Jet Interaction Generator (HIJING) for generating the initial conditions, 
Zhang's Parton Cascade (ZPC) for modeling partonic scatterings, the Lund string fragmentation 
model or a quark coalescence model for hadronization, and A Relativistic Transport (ART) model 
for treating hadronic scatterings, are improved and combined to give a coherent description of the 
dynamics of relativistic heavy ion collisions. We also explain the way parameters in the model are 
determined, and discuss the sensitivity of predicted results to physical input in the model. Com- 
parisons of these results to experimental data, mainly from heavy ion collisions at the Relativistic 
Heavy Ion Collider (RHIC) , are then made in order to extract information on the properties of the 
hot dense matter formed in these collisions. 

PACS numbers: 25.75.-q, 12.38.Mh, 24.10.Lx 



I. INTRODUCTION 



Colliding heavy ions at relativistic energies makes it 
possible to subject nuclear matter to the extreme con- 
dition of large compression, leading to energy densities 
that can exceed that for producing a plasma of decon- 
fined quarks and gluons, that is believed to have existed 
during the first microsecond after the Big Bang. Ex- 
periments at RHIC at the Brookhaven National Labo- 
ratory with center-of-mass energy up to ^/snn = 200 
GeV in Au+Au collisions thus provide the opportunity to 
study the properties of this so-called quark-gluon plasma 
(QGP). At the future Large Hadron Collider (LHC) at 
CERN, which will allow Pb4-Pb collisions at ^/sWn = 5.5 
TeV, the produced quark-gluon plasma will have an even 
higher temperature and a nearly vanishing net baryon 
chemical potential. 

Many observables have been measured at RHIC, such 
as the rapidity distributions of various particles and their 
transverse momentum spectra up to very high transverse 
momentum, the centrality dependence of these observ- 
ables, the elliptic flows of various particles, as well as both 
identical and non- identical two-particle correlations. To 
understand these extensive experimental results, many 
theoretical models have been introduced. It ranges from 
thermal models Q, H H Q based on the assumption 
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of global thermal and chemical equilibrium to hydrody- 
namic models 0, H H H, H 0, El based only on the 
assumption of local thermal equilibrium , an d to trans- 
port models .12, 1 J, Q 111 111 111 El 111 [S [13, [H 

il IM |25|, m that treat non-equilibrium dynamics ex- 
plicitly. The thermal models have been very successful 
in accounting for the yield of various particles and their 
ratios, while the hydrodynamic models are particularly 
useful for understanding the collective behavior of low 
transverse momentum particles such as the elliptic flow 
H H ^1 El ■ Since transport models treat chemical and 
thermal freeze-out dynamically, they are also natural and 
powerful tools for studying the Hanbury-Brown-Twiss in- 
terferometry of hadrons. For hard processes that involve 
large momentum transfer, approaches based on the per- 
turbative quantum chromodynamics (pQCD) using par- 
ton distribution functions in the colliding nuclei have 
been used IH 12^. Also, the classical Yang-Mills the- 
ory has been developed to address the evolution of par- 
ton distribution functions in nuclei at ultra-relativistic 
energies |2l [sl El ^^'^ used to study the hadron rapid- 
ity distribution and its centrality dependence at RHIC 
32, 33, 34]. These problems have also been studied in 
the pQCD based final-state saturation model |35ll36ll37| . 

Although studies based on the pQCD Is^ have shown 
that thermalization could be achieved in collisions of 
very large nuclei and/or at extremely high energy, even 
though the strong coupling constant at the saturation 
scale is asymptotically small, the dense matter created 
in heavy ion collisions at RHIC may, however, not 
achieve full thermal or chemical equilibrium as a re- 
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suit of its finite volume and energy. To address such 
non-equilibrium many-body dynamics, we have devel- 
oped a multi-phase transport (AMPT) model that in- 
cludes both initial partonic and final hadronic interac- 
tions and the transition between these two phases of mat- 
ter ^M,ALMtM,M,M,M £L 48, M, M- The 

AMPT model is constructed to describe nuclear colfisions 
ranging from p + A to A + A systems at center-of-mass 
energies from about ^/SNN = 5 GeV up to 5500 GeV at 
LHC, where strings and minijets dominate the initial en- 
ergy production and effects from final-state interactions 
are important. For the initial conditions, the AMPT 
model uses the hard minijet partons and soft strings from 
the HIJING model. The ZPC parton cascade is then used 
to describe scatterings among partons, which is followed 
by a hadronization process that is based on the Lund 
string fragmentation model or by a quark coalescence 
model. The latter is introduced for an extended AMPT 
model with string melting in which hadrons, that would 
have been produced from string fragmentation, are con- 
verted instead to their valence quarks and antiquarks. 
Scatterings among the resulting hadrons are described 
by the ART hadronic transport model. With parame- 
ters, such as those in the string fragmentation, fixed by 
the experimental data from heavy ion collisions at SPS, 
the AMPT model has been able to describe reasonably 
many of the experimental observations at RHIC. 



In this paper, we give a detailed description of the 
different components of the AMPT model, discuss the 
parameters in the model, show the sensitivity of its re- 
sults to the input to the model, and compare its predic- 
tions with experimental data. The paper is organized 
as follows. In Section^ we describe the different com- 
ponents of the AMPT model: the HIJING model and 
string melting, the ZPC model, the Lund string frag- 
mentation model, and the quark coalescence model used 
for the scenario of string melting, and the extended ART 
model. Tests of the AMPT model against data from pp 
and pp reactions are given in Section llTTl Results from the 
AMPT model for heavy ion collisions at SPS energies are 
discussed in Section llVI for hadron rapidity distributions 
and transverse momentum spectra, baryon stopping, and 
antiproton production. In Section we show results 
at RHIC for hadron rapidity distributions and trans- 
verse momentum spectra, particle ratios, baryon and an- 
tibaryon production, and the production of multistrange 
baryons as well as J/t/j- We further show results from the 
AMPT model with string melting on hadron elliptic flows 
and two-pion interferometry at RHIC. In Section [VII we 
present the predictions from AMPT for hadron rapidity 
and transverse momentum distributions in Pb-|-Pb colli- 
sions at the LHC energy. Discussions on possible future 
improvements of the AMPT model are presented in Sec- 
tion |VH1 and a summary is finally given in Section [VIIII 



II. THE AMPT MODEL 



The AMPT model consists of four main components: 
the initial conditions, partonic interactions, the con- 
version from the partonic to the hadronic matter, and 
hadronic interactions. The initial conditions, which in- 
clude the spatial and momentum distributions of minijet 
partons and soft string excitations, are obtained from 
the HIJING model iHl HI HI iHl- Currently, the 
AMPT model uses the HIJING model version 1.383 [H, 
which does not include baryon junctions '56|. Scatter- 
ings among partons are modeled by Zhang's parton cas- 
cade (ZPC) 18], which at present includes only two-body 
scatterings with cross sections obtained from the pQCD 
with screenin g rn asses. In the default AMPT model 
|39ll4^l4lLl4^l4^l44Ll4Rl47ll^ partons are recombined 
with their parent strings when they stop interacting, and 
the resulting strings are converted to hadrons using the 
Lund string fragmentation model I5a 158, l59l . In the 
AMPT model with string melting [43, |43, llfll , a quark 
coalescence model is used instead to combine partons into 
hadrons. The dynamics of the subsequent hadronic mat- 
ter is described b y a hadronic cascade, which is based on 
the ART model |l4l l25l | and extended to include addi- 
tional reaction channels that are important at high en- 
ergies. These channels include the formation and decay 
of K* resonance and antibaryon resonances, and baryon- 
antibaryon production from mesons and their inverse re- 
actions of annihilation. Final results from the AMPT 
model are obtained after hadronic interactions are ter- 
minated at a cutoff time (tent) when observables un- 
der study are considered to be stable, i.e., when further 
hadronic interactions after tcut will not significantly af- 
fect these observables. We note that two-body partonic 
scatterings at all possible times have been included be- 
cause the algorithm of ZPC, which propagates partons 
directly to the time when the next collision occurs, is 
fundamentally different from the fixed time step method 
used in the ART model. 

In Figs. Hand 12 we show, respectively, the schematic 
structures of the default AMPT model [Hi El [il IH 
m I4II4II4I and the AMPT model with string me\t- 
ing |45lf48l l5(]l| described above. The full source code 
of the AMPT model in the Fortran 77 language and in- 
structions for users are available online at the OSCAR 
website and also at the EPAPS website 61]. The 
default AMPT model is named as version l.a;, and the 
AMPT model with string melting is named as version 
2.y, where value of the integer extension x or y increases 
whenever the source code is modified. Current versions 
of the AMPT model is 1.11 for the default model and 
2.11 for the string melting model, respectively. In the 
following, we explain in detail each of the above four 
components of the AMPT model and the way they are 
combined to describe relativistic heavy ion collisions. 
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Structure of the default AM PT model 
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FIG. 1: (Color online) Illustration of the structure of the 
default AMPT model. 



Structure of AMPT model with string melting ^^^d to fit the same pp and data 
A+B 
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dial density profiles of the two colliding nuclei are taken 
to have Woods-Saxon shapes, and multiple scatterings 
among incoming nucleons are treated in the eikonal for- 
malism. Particle production from two colliding nucleons 
is described in terms of a hard and a soft component. 
The hard component involves processes in which the mo- 
mentum transfer is larger than a cutoff momentum pq 
and is evaluated by the pQCD using the parton distri- 
bution function in a nucleus. These hard processes lead 
to the production of energetic minijet partons and are 
treated via the PYTHIA program. The soft component, 
on the other hand, takes into account non-perturbative 
processes with momentum transfer below po and is mod- 
eled by the formation of strings. The excited strings are 
assumed to decay independently according to the Lund 
JETSET fragmentation model. 

From the pp and pp total cross sections and the ratio of 
CToi/CTtot in the energy range 20 < < 1800 GeV, it has 
been found that the experimental data can be fitted with 
a nucleon-nucleon soft cross section as{s) = 57 mb at 
high energies and po = 2 GeV/c [HJ]. The independence 
of these two parameters on the colliding energy is due to 
the use of the Duke-Owens set 1 for the parton distribu- 
tion function 62] in the nucleon. With different parton 
distribution functions, an energy-dependent pn may be 
needed to fit the same pp and pp data [HE 1^ ■ We note 
that since the number of hard collisions in an A + A colli- 
sion roughly scales as A**/^ and grows fast with colliding 
energy while the number of strings roughly scales as A, 
minijet production becomes more important when the 
energy of heavy ion collisions increases [5ll l65l | . 

Because of nuclear shadowing, both quark and 
giuon [g^ distribution functions in nuclei are different 
from the simple superposition of their distributions in 
a nucleon. This effect has been included in the HIJING 
model via the following impact-parameter-dependent but 
(5^(and flavor)-independent parameterization 52] : 



RA{x,r) = 



l + 1.191n^/^A(x3 - 
1.08(^1/3 



aA{r)- 



1 



ln(A+ 1) 



0.21a;) 

e--'/°-°\(l) 



ART (A Relativistic Transport model for hadrons) 



FIG. 2: (Color online) Illustration of the structure of the 
AMPT model with string melting. 



A. Initial conditions 

1. The default AMPT model 

In the default AMPT model, initial conditions for 
heavy ion collisions at RHIC arc obtained from the HI- 
JING model m m ^ |5J]. In this model, the ra- 



where x is the light-cone momentum fraction of parton a, 
and fa is the parton distribution function. The impact- 
parameter dependence of the nuclear shadowing effect is 
controlled by 



aA{r) = 0A33{A 



1/3 



1)/ 



1 



(2) 



with r denoting the transverse distance of an interacting 
nucleon from the center of the nucleus with radius Ra ~ 
1.2^1/3. Note that there is a modified HIJING model 
which uses a different parameterization for the nuclear 
shadowing that is also flavor-dependent [6^ . 

To take into account the Lorentz boost effect, we have 
introduced a formation time for minijet partons that de- 
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pends on their four momenta |68j |. Specifically, the for- 
mation time for each parton in the default AMPT model 
is taken to have a Lorentzian distribution with a half 
width </ = E/m^, where E and niT are the parton en- 
ergy and transverse mass, respectively. Initial positions 
of formed minijet partons are calculated from those of 
their parent nucleons using straight-line trajectories. 



2. The AMPT model with string melting 

Although the partonic part in the default AMPT 
model includes only minijets from the HIJING model, its 
energy density can be very high in heavy ion collisions 
at RHIC. As shown in Fig. O for the time evolutions of 
the energy and number densities of partons and hadrons 
in the central cell of central (6 = fm) Au+Au colli- 
sions at y/SNN — 200 GeV in the center-of-mass frame, 
the partonic energy density during the first few fm/c's of 
the collision is more than an order of magnitude higher 
than the critical energy density 1 GeV/fm^) for the 
QCD phase transition, similar to that predicted by the 
high density QCD approach [6^. The sharp increase in 
energy and number densities at about 3 fm/c is due to 
the exclusion of energies that are associated with the ex- 
cited strings in the partonic stage. Keeping strings in 
the high energy density region |7(]| thus underestimates 
the partonic effect in these collisions. We note that the 
central cell in the above calculation is chosen to have 
a transverse radius of 1 fm and a longitudinal dimension 
between —0.5t and 0.5t, where time t starts when the two 
nuclei are fully overlapped in the longitudinal direction. 

To model the above effect in high energy density re- 
gions, we extend the AMPT model to include the string 
melting mechanism [s^, i.e, all excited strings, 

that are not projectile and target nucleons without any 
interactions, are converted to partons according to the 
flavor and spin structures of their valence quarks. In 
particular, a meson is converted to a quark and an anti- 
quark, while a baryon is first converted to a quark and 
a diquark with weights according to relations from the 
SU(6) quark model [71|, and the diquark is then decom- 
posed into two quarks. The quark and diquark masses 
are taken to be the same as in the PYTHIA program [s^ , 
e.g. rriu = 5.6 MeV/c^, rud = 9.9 MeV/c^ and = 199 
MeV/c^. We further assume that the above two-body 
decomposition is isotropic in the rest frame of the par- 
ent hadron or diquark, and the resulting partons do not 
undergo scatterings until after a formation time given by 
tf = Eh /rn^ jj, with Eh and itit.h denoting the energy 
and transverse mass of the parent hadron. Similar to 
minijet partons in the default AMPT model, initial po- 
sitions of the partons from melted strings are calculated 
from those of their parent hadrons using straight- line tra- 
jectories. 

The above formation time for partons is introduced 
to represent the time needed for their production from 
strong color fields. Although we consider hadrons be- 
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FIG. 3: Energy and number densities of minijet partons and 
formed hadrons in the central cell as functions of time for 
central (b=0 fm) Au-|-Au collisions at y^ijvjv = 200 GeV from 
the default AMPT model, where the energy stored in the 
excited strings is absent in the parton stage and is released 
only when hadrons are formed. 



fore string melting as a convenient step in modeling the 
string melting process, choosing a formation time that 
depends on the momentum of the parent hadron ensures 
that partons from the melting of same hadron would have 
the same formation time. The advantage of this choice 
is that the AMPT model with string melting reduces to 
HIJING results in the absence of partonic and hadronic 
interactions as these partons would then find each other 
as closest partners at the same freeze-out time and thus 
coalesce back to the original hadron. We note that the 
typical string fragmentation time of about 1 fm/c is not 
applied to the melting of strings as the fragmentation 
process involved here is considered just as an intermedi- 
ate step in modeling parton production from the energy 
field of the strings in an environment of high energy den- 
sity. 



B. Parton cascade 

In the transport approach, interactions among partons 
are described by equations of motion for their Wigner dis- 
tribution functions, which describe semi-classically their 
density distributions in phase space. These equations 
can be approximately written as the following Boltzmann 
equations: 



p'^9^/a(x,p,t) 



5 



E E n 



m bi,b2,---,b„ 



-/fc.(x, Pbi,0 



E E /n7;T^iM„_p 



n ci ,C2 , ■ ■ -.Ct, 



1 (27r)32ii;,^ 



\k=l 



1 = 1 



E '^'^^'^^(P - Pbq) + E '^"^.'^^(P - PcJ 
9—1 r— 1 



(3) 



In the above, /q(x. p,i) is the distribution function of 
parton type a at time t in the phase space, and Af„j^„ 
denotes the matrix element of the muhi-parton interac- 
tion m —I- n. If one considers only two-body interactions, 
these equations reduce to 



P^df,f{x,p,t) (X J (T/(xi,Pl,i)/(x2,P2,t), 



(4) 



where a is the cross section for partonic two-body scat- 
tering, and the integral is evaluated over the momenta of 
other three partons with the integrand containing factors 
such as a delta function for momentum conservation. 

The Boltzmann equations are solved using Zhang's 
parton cascade (ZPC) in which two partons undergo 
scattering whenever they approach each other with a 
closest distance smaller than ^JuJt^. At present, ZPC in- 
cludes only parton two-body scattering such as gg — > gg 
with cross sections calculated from the pQCD. For gluon 
elastic scattering, the leading-order QCD gives 



dt 
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where Us is the strong coupling constant, and s, t and u 
are standard Mandelstam variables for elastic scattering 
of two partons. The second line in the above equation 
is obtained by keeping only the leading divergent terms. 
Since the scattering angle ranges from to 7r/2 for iden- 
tical particles, one then has 



doa 



dt 



9TTai 



(6) 



if the scattering angle is between and tt. 

The singularity in the total cross section can be regu- 
lated by a Debye screening mass /i, leading to 



dfj, 



gg 



dt 



2(i-/i2)2' 

97rQ;2 1 

2^2 l+fj,2/s' 



(7) 



The screening mass fj, is generated by medium effects and 
is thus related to the parton phase-space density. For the 



partonic system expected to be formed in Au+Au colli- 
sions at RHIC, the value of fj, is on the order of one inverse 
fermi 18]. For massless partons in a plasma at temper- 
ature T, their average colliding energy is y/s ^ vTST, 
thus n < y/s for ^ = 3 fm~^ leads to the requirement 
T > 141 MeV. Since s > fi^ generally holds in hot QGP, 
the following simplified relation between the total parton 
elastic scattering cross section and the screening mass is 
used in the ZPC Fl 



'gg 



9'Kai 



(8) 



A value of 3 fm^^ for the screening mass fi thus leads to a 
total cross section of about 3 mb for the elastic scattering 
between two gluons. By changing the value of the screen- 
ing mass fl, different cross sections can be obtained, and 
this will be used in studying the effect of parton cross sec- 
tions in heavy ion collisions at RHIC. This cross section 
is used in AMPT not only in the default model, which 
includes only scatterings of minijet gluons, but also in 
the string melting model, which only includes scatter- 
ings of quarks/antiquarks of all flavors. We have there- 
fore neglected in the latter case the difference between 
the Casimir factors for quarks and gluons. 

We note that minijet partons produced from hard scat- 
terings in the HIJING model can lose energy by gluon 
splitting and transfer their energies to nearby soft strings. 
In the AMPT model, this so-called jet quenching in the 
HIJING model is replaced by parton scatterings in ZPC. 
Since only two-body scatterings are included in ZPC, 
higher-order contributions to the jet energy loss are still 
missing in the AMPT model. 



C. Hadronization 

Two different hadronization mechanisms are used in 
the AMPT model for the two different initial conditions 
introduced in Sec. Ill Al In the default AMPT model, 
minijets coexist with the remaining part of their parent 
nucleons, and they together form new excited strings af- 
ter partonic interactions. Hadronization of these strings 
are described by the Lund string model. In the AMPT 
model with string melting, these strings are converted to 
soft partons, and their hadronization is based on a simple 
quark coalescence model, similar to that in the ALCOR 
model |73j. 



1. Lund string fragmentation for the default AMPT model 

Hadron production from the minijet partons and soft 
strings in the default AMPT model is modeled as follows. 
After minijet partons stop interacting, i.e., after they no 
longer scatter with other partons, they are combined with 
their parent strings to form excited strings, which are 
then converted to hadrons according to the Lund string 
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fragmentation model |53, |53| . In the Lund model as im- 
plemented in the JETSET/PYTHIA routine 59], one as- 
sumes that a string fragments into quark-antiquark pairs 
with a Gaussian distribution in transverse momentum. 
A suppression factor of 0.30 is further introduced for the 
production of strange quark-antiquark pairs relative to 
that of light quark-antiquark pairs. Hadrons are formed 
from these quarks and antiquarks by using a symmetric 
fragmentation function 57, 58] . Specifically, the trans- 
verse momentum of a hadron is given by those of its 
constituent quarks, while its longitudinal momentum is 
determined by the Lund symmetric fragmentation func- 
tion jlll 

/(z) cx z-i(l - z)"exp(-6 TO^/z), (9) 

with z denoting the light-cone momentum fraction of the 
produced hadron with respect to that of the fragment- 
ing string. The average squared transverse momentum is 
then given by 



J p\f{z)(fpj_dz 

J f{z)cPp^dz 
jj;max _ 2:)aexp(— 6 m? / z)dz 
b /g^"""^ (1 — z)" exp(— 6 m? / z)dz 
For massless particles, it reduces to 



(pi) 



1 z(l - zYdz _ 1 

b j'l{l-zYdz ^ b{2 + a)' 



(10) 



(11) 



Since quark-antiquark pair production from string 
fragmentation in the Lund model is based on the 
Schwinger mechanism ff^ for particle production in 
strong field, its production probability is proportional 
to exp(— tttti^/k), where k is the string tension, i.e., 
the energy in a unit length of string. Due to its large 
mass, strange quark production is suppressed by the 
factor e"'^'-™'"™"^/", compared to that of light quarks. 
Also, the average squared transverse momentum of pro- 
duced particles is proportional to the string tension, i.e., 
(p^) oc K. Comparing this with Eq. one finds that 

the two parameters a and b in the Lund fragmentation 
function are approximately related to the string tension 

by 

1 , X 

i^oc— -. 12) 

6(2 + a) ^ ' 

After production from string fragmentation, hadrons 
are given an additional proper formation time of 0.7 fm/c 
fr^ . Positions of formed hadrons are then calculated 
from those of their parent strings by following straight- 
line trajectories. 



2. Quark coalescence for the AMPT model with string 
melting 

After partons in the string melting scenario stop inter- 
acting, we model their hadronization via a simple quark 



coalescence model by combining two nearest partons into 
a meson and three nearest quarks (antiquarks) into a 
baryon (antibaryon) . Since the invariant mass of com- 
bined partons forms a continuous spectrum instead of 
a discrete one, it is generally impossible to conserve 4- 
momentum when partons are coalesced into a hadron. 
At present, we choose to conserve the three-momentum 
during coalescence and determine the hadron species ac- 
cording to the flavor and invariant mass of coalescing 
partons [t^I- For pseudo-scalar and vector mesons with 
same flavor composition, the meson with mass closer to 
the invariant mass of coalescing quark and antiquark pair 
is formed. E.g., whether a 7r~ or a p~ is formed from the 
coalescence of a pair of u and d quarks depends on if the 
invariance mass of the quarks is closer to the tt" mass or 
the centroid of p mass. The same criterion applies to the 
formation of octet and decuplet baryons that have same 
flavor composition. It is more complicated to treat the 
formation probabilities of flavor-diagonal mesons such as 
tt'^ and rj in the pseudo-scalar meson octet, and and 
w in the vector meson octet. Neglecting the mixing of rj 
meson with the ss state, we take the following approach 
for these flavor-diagonal mesons within the SU(2) flavor 
space. For tt^ formation from a uu or dd pair, the proba- 
bility PttO is determined from the average of the numbers 
of formed 7r+ and 7r~ mesons divided by the total number 
of uu and dd pairs. Thus the total number of tt*^, n^o, is 
determined by applying the probability P^o to each uu or 
dd pair. The probability PpO for forming a meson from 
a uu or dd pair is determined by a similar procedure. Af- 
ter sorting all uu or dd pairs according to their invariant 
masses, the lightest Ut^o pairs are assigned as 7r° mesons. 
The rest of uu or dd pairs form mesons according to 
the probability of PpO /(I — Ptt"), and the remaining pairs 
form Lo and 77 mesons with equal probabilities. 

The above quark coalescence model includes the for- 
mation of all mesons and baryons listed in the HIJING 
program except 77', S* and which are not present 
in our hadronic transport model as well as Kg and 
states. The resulting hadrons are given an additional for- 
mation time of 0.7 fm/c in their rest frame before they are 
allowed to scatter with other hadrons during the hadron 
cascade. As partons freeze out dynamically at different 
times in the parton cascade, hadron formation from their 
coalescence thus occurs at different times, leading to the 
appearance of a coexisting phase of partons and hadrons 
during hadronization. 



D. Hadron cascade 

In the AMPT model, the following hadrons with all 
possible charges are explicitly included: tt, p, uj, rj, K, 
K*, and (j) for mesons; iV, A, 7V*(1440), iV*(1535), A, 
S, S, and VI for baryons and corresponding antibaryons. 
Many other higher resonances are taken into account 
implicitly as intermediate states in scatterings between 
the above particles 0, |2^. Interactions among these 
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hadrons and corresponding inverse reactions are included 
as discussed in following subsections. 



1. The ART model 

Hadron cascade in the AMPT model is based on the 
ART model IZ^i which is a relativistic transport 
model originally developed for heavy ion collisions at 
AGS energies. The ART model includes baryon-baryon, 
baryon-meson, and meson-meson elastic and inelastic 
scatterings. It treats explicitly the isospin degrees of 
freedom for most particle species and their interactions, 
making it suitable for studying isospin effects in heavy 
ion collisions ,79i] . Since it includes mean-field poten- 
tials for nucleons and kaons, the ART model can also be 
used for studying the effect due to the hadronic equation 
of state. Resonances such as p and A are formed from 
pion-pion and pion-nucleon scattering, respectively, with 
cross sections given by the standard Breit-Wigner form, 
and they also decay according to their respective widths. 
In all calculations presented in this study, the masses 
and widths of resonances are taken to be their values in 
the vacuum, i.e., effects due to possible modifications in 
dense hadronic matter are neglected. Also, we have 
turned off the potentials in the AMPT model as their 
effects are much less important than scatterings in high 
energy heavy ion collisions such as at SPS and RHIC. 

For baryon-baryon scatterings, the ART model in- 
cludes the following inelastic channels: NN ^ N{AN*), 
NN ^ A(AiV*(1440)), NN ^ NNi-npuj), {NA)A ^ 
NN*, and A7V*(1440) ^ A^A^*(1535). In the above, 
N* denotes either iV*(1440) or iV*(1535), and the sym- 
bol (AiV*) denotes a A or an N*, Also included are 
reaction channels relevant for kaon production, i.e., 
{NAN*){NAN*) {NA){AY.)K. Details on their 
cross sections and the momentum dependence of reso- 
nance widths can be found in the original ART model 

For meson-baryon scatterings, the ART model includes 
the following reaction channels for the formation and de- 
cay of resonances: TriV ^ {AN*{U40) 7V*(1535)), and 
r]N <-> A^*(1535). There are also elastic scatterings such 
as {ttp){NAN*) {ttp){NAN*). As an example, the 
cross section for the elastic scattering of tt'^N is evalu- 
ated by including heavier baryon resonances with masses 
up to 2.0 GeV/c^ as intermediate states using the Breit- 
Wigner form but neglecting interferences between the 
amplitudes from different resonances [13. The ART 
model further includes inelastic reaction channels such 
as TriV <-> {-KprfjA and kaon production channels such 
as {npLorj){N AN*) <-> K{KYi). Kaon elastic scatterings 
with nucleons and baryon resonances are included with 
a constant cross section of 10 mb Antikaon elastic 
scatterings with nucleons and inelastic channels, such as 
K{NAN*) ^ 7r(Ai;), are included using parame- 
terized experimental data [82j. Also included are kaon 
production channels involving three-body final states. 



{TTpuj){NAN*) KKN |8l|. Because of the difficulty 
associated with the three-body kinematics, the inverse 
kaon annihilation reactions of the above channels are ne- 
glected. 

For meson-meson interactions, the ART model in- 
cludes both elastic and inelastic tttt interactions, with 
the elastic cross section consisting of p meson forma- 
tion and the remaining part treated as elastic scatter- 
ing. Kaon production from inelastic scatterings of light 
mesons is included via the the reactions (ttt^) (ttt^) <-> KK 
and {puj)(puj) <-!■ KK. Kaon or antikaon elastic scatter- 
ings with mesons in the SU(2) multiplets except the pion 
are included using a constant cross section of 10 mb |14| . 
while the kaon-pion elastic scattering is modeled through 
the K* resonance 



2. Explicit inclusion of K* mesons 

The original ART model Q includes the K* reso- 
nance implicitly through elastic ttK scattering with the 
standard Breit-Wigner form for the cross section [83l| . 
i.e., cr^K = 60 mb/ [l + - mK- f /Vl^,]. Since the 
K* meson not only enhances elastic scattering between 
pion and kaon but also adds to strange particle produc- 
tion through its addition to the strangeness degeneracy, 
which becomes important when the hadronic matter is 
highly excited, K* and K* are included explicitly in the 
hadronic phase of the AMPT model addition 
to its formation from t:K scattering and its decay, elas- 
tic scatterings of K* with {poJVi) are included using same 
constant cross section of 10 mb as those for kaons. In- 
elastic reaction channels of {'Krj){pijj) <-> K*K or K*K 
and ttK ^ K*{pu) are also included ( 84] . 



3. Baryon-antibaryon annihilation and production 

In heavy ion collisions at or above SPS energies, an- 
tibaryon production becomes significant and needs to be 
treated explicitly during hadron cascade. The AMPT 
model initially only includes NN annihilation It 
was later extended to include {NAN*){NAN*) annihi- 
lation and also the inverse reactions of baryon-antibaryon 
pair production from mesons _42J. 

The total cross section for pp annihilation is known 
empirically, and the data has been parameterized as [ssj ] 



app = 67 mh/p^J^, 



(13) 



where piab in GeV/c is the proton momentum in the rest 
frame of the antiproton. Following Ref. [l^ . a maximum 
cross section of 400 mb is imposed at low piab- Using 
phase-space considerations jSg, the branching ratios of 
pp annihilation to different multi-pion states are deter- 
mined according to |86j | 



M„(Vi) ^ C 



1 

6^ 



(4n-4)!(2n- 1) 
{2n- l)!2(3n- 4)! 



,(14) 
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where ^/s is the center-of-mass energy of the proton and 
antiproton, and n is the number of pions in the final 
state. The dominant final states at moderate energies 
then involve several pions. For example, the branch- 
ing ratios at piab = 4 GeV/c are 0.033, 0.161, 0.306, 
0.286, 0.151, 0.049 and 0.011, respectively, for n from 3 
to 9. For baryon-antibaryon annihilation channels involv- 
ing baryon resonances A or iV*, their annihilation cross 
sections and branching ratios are taken to be the same 
as for pp annihilation at same center-of-mass energy. 

To include the inverse reactions that produce baryon- 
antibaryon pairs during hadron cascade, which currently 
only treats scattering of two particles, we have further 
assumed that the final state of three pions is equivalent 
to a TTyO state; the four-pion final state is equivalent to 
pp and TTw with equal probabilities; the five-pion state is 
equivalent to puj] and the six-pion state is equivalent to 
Ljuj. The cross sections for baryon-antibaryon pair pro- 
duction from two mesons are then obtained from detailed 
balance relations. As shown later in the paper fFig. I18|l . 
the above approximate treatment of antibaryon annihila- 
tion and production via two-meson states gives a satisfac- 
tory description of measured antiproton yield in central 
Pb+Pb colhsions at SPS. 



^. Multistrange baryon production from 
strangeness-exchange reactions 



strangeness-exchange reactions such as 
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FIG. 4: (Color online) Isospin-averaged cross sections for 
(multi)strange baryon production as functions of center-of- 
mass energy of interacting antikaon and nucleon/hyperon. 



Productions of multistrange baryons such as :^ and 
are included in the AMPT model Q through the 



i^(AE) 



ttS, KE 



(15) 



Since there is no experimental information on their cross 
sections, we have assumed that the matrix elements for 
K{AI]) ttS and KE ttQ are the same as that for the 
reaction KN — > ttS [s^I at the same amount of energy 
above corresponding thresholds. The isospin-averaged 
cross section for the reaction KN — > ttE can be related 
to the cross sections for the reactions K~p — + E^tt" and 
K^n — > E'^TT^, which are known empirically and have 
been parameterized in Ref. [s^l, by 



'KN- 



;(fi<'-p^E"7rO + fK-n-fSOu--)- (16) 



In Fig.^ this cross section as well as the isospin-averaged 
cross sections for other multistrange baryon production 
reactions are shown as functions of the center-of-mass 
energy above the threshold values of the interacting an- 
tikaon and nucleon/hyperon. The cross sections for 
the inverse multistrange baryon destruction reactions 
are then determined by detailed balance relations. We 
note that these cross sections are comparable to those 
predicted by the coupled-channel calculations based on 
the SU(3) invariant hadronic Lagrangian with empirical 
masses and coupling constants [Sg. With these cross 
sections, the ART transport model is able to describe 
the measured S production in heavy ion collisions at the 
AGS energies jS^- We note that, for strange baryons A, 
S, O and their anti-particles, only their interactions with 
mesons have been included, while their annihilations by 
baryons have not been included in the AMPT model at 
present. 



5. (p meson production and scattering 

The AMPT model also includes (j> meson formation 
from and decay to kaon-antikaon pair with the forma- 
tion cross section given by the standard Breit-Wigner 
form Inelastic scatterings of the phi meson include 
baryon-baryon channels, {NAN*){NAN*) -> (j)NN, 
and meson-baryon channels, {Trp){NAN*) ^ 0(A^AiV*), 
where the cross sections for the forward-going reactions 
are taken from the one-boson-exchange model [o^l . The 
meson-baryon channels also include K{AY,) (j)N with 
cross section taken from a kaon-exchange model [9ll |. 

Phi meson scatterings with mesons such as tt, p, /iT and 
(f) have been studied before, and the total coUisional width 
was found to be less than 35 MeV/c^ 92] . A recent 
calculation based on the Hidden Local Symmetry La- 
grangian '9^ shows, however, that the coUisional rates 
of (j) with pseudo-scalar (tt, K) and vector (p, w, if*, 
0) mesons are appreciably larger. Assuming that the 
matrix elements are independent of center-of-mass en- 
ergy, we have included all these possible reactions, i.e., 
(/)(7rpw) ^ {KK*){KK*), and (t){KK*) <-> {■Kpw){KK*), 
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with cross sections determined from the partial collisional 
widths given in Ref. ^9§\. The cross section for the elas- 
tic scattering of the (f) meson with a nucleon is set to 8 
mb while the (j> meson elastic cross section with a meson 
is set to 5 mb. The value of 8 mb is the (l)N total cross 
section |91l estimated from the (p meson photoproduction 
data jOJI and thus represents the upper bound on the (fi 
meson elastic cross section with a nucleon; quark count- 
ing then gives 5 mb as the upper bound on the (f) meson 
elastic cross section with a meson. Note that, in most 
calculations of our previous study on (j) meson produc- 
tions at SPS and RHIC the elastic cross section for 
(j) meson scattering with a nucleon was taken to be 0.56 
mb as extracted in Ref. ||95|| using the vector meson dom- 
inance model and the older (j) meson photoproduction 
data, while based on results of Ref. [93| the (t>K elastic 
cross section was extracted to be about 2 mb !47| , which 
was then used as the (f) meson elastic cross section with 
other mesons [96|. 

6. Other extensions 

Other extensions of the ART model have also been 
made in the hadronic phase of the AMPT model. An- 
tibaryon resonances such as A and N* have been included 
explicitly with their formations, decays, and scatterings 
analogous to those of baryon resonances • Also, in- 

elastic meson-meson collisions such as tttt ^ pp have been 
added, and elastic scatterings between tt and {pwr]) have 
been included with cross sections taken to be 20 mb. To 
address chemical equilibration of rj mesons, which affects 
the height (or the A parameter) of the correlation func- 
tions in two-pion interferometry, inelastic scatterings of 
77 meson with other mesons have also been included with 
a constant cross section of 5 mb which is roughly 
in line with recent theoretical predictions based on the 
Hidden Local Symmetry Lagrangian 

III. RESULTS FOR pp AND pp COLLISIONS 

The default AMPT results with no popcorn mechanism 
(see discussions on the popcorn mechanism in Sec lIVB)) 
for pp and pp collisions are essentially the same as the 
results from the HIJING model. In this section, we com- 
pare the results from the default AMPT model with or 
without the popcorn mechanism against available data 
from pp and pp collisions , where HIJING values for 
the a and b parameters, a = 0.5 and b = 0.9 GeV~^, are 
used in Eq. ||^ for the Lund fragmentation function. 

A. Rapidity distributions 

In Fig. [SI the AMPT results on charged particle pseu- 
dorapidity distribution are compared with the UA5 data 
for pp colhsions at Vs = 200 GeV |9g|. Since the AMPT 



results given by the solid or the long-dashed curves rep- 
resent all inelastic events with no trigger conditions, 
they should be compared with the UA5 inelastic data. 
The AMPT non-single-diffractive (NSD) results have in- 
cluded the UA5 NSD trigger by requiring events to have 
at least one charged particle each in both ends of the 
pseudorapidity intervals 2 < jT^j < 5.6 We see that 
the AMPT model with or without the popcorn mecha- 
nism agrees reasonably with both NSD and inelastic data. 
The kaon rapidity distribution from the AMPT model is 
compared with the UA5 NSD data in Fig. El and both re- 
sults with or without the popcorn mechanism also agree 
reasonably with the data. 
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FIG. 5: (Color online) Pseudorapidity distributions of 
charged particles for pp collisions at ^/s = 200 GeV. 

For pp collisions at Piab — 400 GeV/c, the rapidity dis- 
tributions of pions, kaons, protons and antiprotons are 
shown, respectively, in Fig.[7| Fig.|Sl Fig.^ and Fig. 1101 
Curves with circles are me asur ed cross sections from the 
LEBC-EHS Collaboration |lO0l |. while for comparison the 
AMPT results have been scaled up by the inelastic cross 
section at this energy {cFinei — 32 mb). We see that simi- 
lar descriptions of the charged particle data are obtained 
with and without the popcorn mechanism. Including the 
popcorn mechanism gives, however, a better description 
of measured pion, kaon and antiproton yields and the 
shape of the proton rapidity distribution. We have thus 
included the popcorn mechanism in all AMPT calcula- 
tions. 



B. Transverse momentum spectra 

The transverse momentum spectra of charged pions 
and protons for pp collisions at Piab = 24 GeV/c from 
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FIG. 6: (Color online) Same as Fig. I^for kaons. 

30 I > ^ > 




FIG. 7: (Color online) Rapidity distributions of pions for pp 
collisions at Piab = 400 GeV/ c. Circles are data from the 
LEBC-EHS Collaboration [Tooll. 

the AMPT model are shown in Fig. ^] They are seen 
to re prod uce reasonably well the experimental data from 
Ref. [loH . 

At the Tevatron energy of y/s — 1.8 TeV, results from 
the AMPT model for the transverse momentum spectra 
of pions, kaons and antiprotons are shown in Fig.|12land 
compared with data from the E735 Collaboration [102ij . 
The AMPT NSD results have included the trigger by re- 
quiring events to have at least one charged particle each 
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FIG. 8: (Color online) Same as Fig.|7|for kaons. 
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FIG. 9; (Color online) Same as Fig.|7|for protons. 



in b oth ends of the pseudorapidity intervals 3 < \ri\ < 4.5 
1 103']. We see that measured momentum spectra except 
for antiprotons are reproduced reasonably well by the 
AMPT model. Note that the E735 pr spectrum data 
shown in Fig. 1121 have been averaged over rapidity y from 
weighing each track by the rapidity range of the spec- 
trometer |l02j . thus they are for cP N / {T.-KpTdpTdy) in- 
stead for (P N / {2'KpT^dp^dr]) even though the E735 spec- 
trometer covers the acceptance of —0.36 < rj < 1.0. 
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FIG. 10: (Color online) Same as Fig. |7|for antiprotons. 




^0 0.5 1 1.5 

FIG. 11: (Color online) Transverse momentum spectra of 
charged pions and protons for pp collisions at Piab ~ 24 GeV/c 
with data from Ref. ilOlj . 



FIG. 12: (Color online) Transverse momentum spectra of 
charged pions, charged kaons and antiprotons from pp col- 
lisions at ^ = 1.8 TeV. 



below y/s — 1 TeV have included the UA5 NSD trigger 
by requiring events to have at least one charged particle 
each in both ends of the pseudorapidity intervals of 2 < 
|?7| < 5.6 [93. The agreement with the Tevatron data is 
reasonable. 
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C. Energy dependence 

Fig. [T^ shows the energy dependence of dNch/drj at 
r] = for pp and pp collisions together with UA5 and 
CDF data at the Tevatron 99, 104]. The AMPT NSD 
results above = 1 TeV have included the CDF NSD 
trigger by requiring events to have at least one charged 
particle each in both ends of the pseudorapidity intervals 
of 3.2 < \t]\ < 5.9 [l03, and the AMPT NSD results 



10 1,0 10 

Vs (GeV) 

FIG. 13: Energy dependence of dNch/drj at mid- 
pseudorapidity for pp and pp collisions. 

The energy dependence of the full phase space A'+/7r+ 
ratio for pp collisions is shown in Fig. E| together with 
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the data compiled in Fig. 7 of Ref. [105| . It is seen that 
the AMPT model reproduces the data reasonably well. 
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FIG. 14: Energy dependence of the /tv'^ ratio in full phase 
space for pp collisions. 

The energy dependence of the mean transverse mo- 
menta of pions, kaons and antiprotons are shown in 
Fig. ^1 where filled circl es re present the NSD data from 
the E735 Collaboration jlOfi| | and open circles represent 
data from the CERN Intersecting Storage Rings (ISR) 
[ToTj . Note that the AMPT NSD results shown in Fis.HB 
have included the E735 NSD trigger |l03l | by requiring 
events to have at least one charged particle each in both 
ends of the pseudorapidity intervals 3 < Ir]] < 4.5. How- 
ever, the trigger condition for the CERN ISR data is 
different. From the comparison with the AMPT inelas- 
tic results (solid curves), we see that NSD triggers have 
larger effects at lower energies. 



IV. RESULTS AT SPS ENERGIES 

To make predictions for heavy ion collisions at RHIC, 
we first use the AMPT model to study heavy ion col- 
lisions at SPS. In particular, parameters in the AMPT 
model are determined by fitting the experimental data 
from central Pb-|-Pb collisions at the laboratory energy 
of 158^ GeV, corresponding to a center-of-mass energy 
of about y/sNN = 17.3 GeV. 
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FIG. 15: (Color online) Energy-dependence of the mean 
transverse momenta of pions, kaons and antiprotons for pp 
and pp collisions. 
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FIG. 16: (Color online) Rapidity distributions in central {b < 
3 fm) Pb-I-Pb collisions at ^sjvjv = 17.3 GeV. Circles and 
squares are experimental data, while dashed curves are results 
from the AMPT model using default a and b parameters as 
in the HIJING model. 



A. Rapidity distributions 

In Fig. El we show the rapidity distributions of neg- 
atively charged particles (upper left panel), net-protons 
and antiprotons (upper right panel) , charged pions (lower 



left panel) , and charged kaons (lower right panel) in cen- 
tral {b < 3 fm) Pb+Pb collisions at SPS, obtained from 
the AMPT model with the default values of a = 0.5 and 
b = 0.9 GeV-2 in the HIJING model for the string frag- 
mentation function. Compared with experimental data 
for 5% most central Pb+Pb collisions from the NA49 
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Collaboration [lOa |108|, ll09| , it is seen that the AMPT 
model with these a and h parameters under-predicts both 
the negatively charged particle 39] and kaon multiplici- 
ties. 
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FIG. 18: (Color online) Same as Fig. ^3 but with dashed 
curves for results from the AMPT model without partonic 
and hadronic final-state interactions. 



FIG. 17: (Color online) Same as Fig. I16l but with solid curves 
for results from the AMPT model using modified a and b 
parameters in the Lund string fragmentation function. 

To increase the particle multiplicity in the AMPT 
model, we vary the a and b parameters, and find that 
the experimental data can be reasonably reproduced with 
a = 2.2 and b — 0.5 GeV~^ as shown in Fig. El As seen 
from Eq. @, a larger value of a corresponds to a softer 
fragmentation function, i.e., it leads to a smaller average 
transverse momentum for produced hadrons and thus in- 
creases the particle multiplicity. The modified parame- 
ters also affect the string tension. According to Eq. (|12|l , 
it is now 7% larger, and this leads to a strangeness sup- 
pression factor of 0.33 instead of the default value of 0.30. 
Since the HIJING model with default a and b parame- 
ters reproduces the charged particle multiplicities in pp 
and pp collisions, the AMPT model with these parame- 
ters can probably describe peripheral heavy ion collisions. 
The a and b values in the AMPT model is thus expected 
to depend on the centrality of heavy ion collisions, but 
this has not been studied. 

Since the probability for minijet production in the HI- 
JING model is very small in collisions at SPS energies 
(only about 4 minijet partons for a central Pb+Pb event 
on average), the partonic stage in the default AMPT 
model does not play any role for most observables in these 
collisions. Final-state hadronic scatterings are, however, 
important, and their effects are illustrated in Fig.^J It 
is seen that final-state interactions reduce the pion and 
antiproton yields, but increase strangeness production, 
e.g., the sum of and A as well as the sum of if + and 



A are both increased by about 20%. Overall, the kaon 
yield in the default AMPT model is larger than that in 
the HIJING model by about 40% due to combined effects 
of modified Lund string fragmentation and final-state in- 
teractions. 



B. Baryon stopping and antiproton production 

Understanding baryon stopping in relativistic heavy 
ion collisions is important as it is related to the total en- 
ergy deposited in the produced hot dense matter during 
the collisions. The observed relatively large baryon stop- 
ping in these collisions has led to the novel suggestion of 
gluon junction transport in an initial excited baryon and 
its subsequent decay into a slowly moving baryon and 
three leading beam mesons jSfij. In the AMPT model, 
we have taken instead a more phenomenological approach 
to baryon stopping in relativistic heavy ion collisions by 
including the popcorn mechanism for baryon-antibaryon 
pair production from string fragmentation. The popcorn 
mechanism introduces two additional baryon production 
channels, i.e., the BB and BMB configurations in the 
Lund fragmentation model, which are controlled by two 
parameters in the JETSET/PYTHIA routine 59] used 
in the HIJING model. The first parameter MSTJ(12) is 
changed from 1 as set in the default HIJING model to 3 in 
the AMPT model 40, 42] in order to activate the popcorn 
mechanism, and the second parameter PARJ(5) controls 
the relative percentage of the BB and BMB channels. 
We find that with equal probabilities for the BB and 
BMB configurations the net-baryon rapidity distribu- 
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tion at SPS can be reproduced as shown in Fig. 1171 [i^. 
Without the popcorn mechanism, as in the default HI- 
JING model, the net-baryon rapidity distribution would 
peak at a larger rapidity • We also see from Fig. ITI^ 
that the antiproton yield at SPS is sensitive to the an- 
tibaryon annihilation and production channels discussed 
in Sec. Ill D 31 Without these reactions in the hadronic 
phase, the antiproton yield is too high compared to pre- 
liminary NA49 data. 



C. Transverse momentum spectra 
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FIG. 19: (Color online) Transverse momentum spectra at 
mid-rapidity from the AMPT model with and without includ- 
ing hadron cascade compared with experimental data from 
the NA44 collaboration for central Pb-|-Pb collisions at SPS. 



For the transverse momentum spectra, results from 
the default AMPT model are shown by solid curves in 
Fig. ll9l for mid-rapidity pions, kaons, and protons in cen- 
tral {b < 3 fm) Pb-I-Pb collisions at the SPS energy of 
y/sNN — 17.3 GeV. Compared with experimental data 
from NA44, given by solid diamonds, the AMPT model 
gives a reasonable description up to transverse mass of 
about 1 GeV/c^ above the particle mass. Without in- 
cluding rescatterings in the hadronic phase, the inverse 
slope parameters for the transverse mass spectra of kaons 
and protons from the AMPT model, given by dashed 
curves, are significantly reduced due to the absence of 
transverse flow that is induced by final-state interactions. 



V. RESULTS AT RHIC ENERGIES 



A. Rapidity distributions 




FIG. 20; (Color online) Rapidity distributions for Au-f-Au col- 
lisions at ^/sWn = 130 and 200 GeV. Circles are the PHOBOS 
data for 6% most central collisions or the BRAHMS data for 
5% most central collisions at y^ijvjv = 200 GeV, and curves 
are results from the default AMPT model for 6 < 3 fm. 

Since the number of strings associated with soft in- 
teractions in the HIJING model depends weakly on the 
colliding energy, and the atomic number of Au is close 
to that of Pb, we use the same modified parameters in 
the Lund string fragmentation model used for Pb-fPb 
collisions at SPS energies for Au+Au collisions at RHIC 
energies. In Fig. 1201 results for central (6 < 3 fm) Au+Au 
collisions at center-of-mass energies of y^sjvw = 130 GeV 
(dashed curves) and 200 GeV (solid curves) are shown 
together with data from the PHOBOS Collaboration 
ESEiiiiilll and the BRAHMS Collaboration 11^. 
We find that measured total charged particle pseudo- 
rapidity distributions at both energies are roughly re- 
produced. More detailed comparisons on pseudorapidity 
distributions at different centralities and different RHIC 
energies have bee n carried out by the BRAHMS Collab- 
oration where results from the default AMPT 
model are compatible with the data. However, compared 
to central (top 5%) BRAHMS data, the AMPT model 
tends to over-predict the height of the rapidity distribu- 
tions of charged pions, kaons, protons and antiprotons. 
Note that the BRAHMS data on proton and antiprotons 
at y = in Fig. [201 have been corrected for feed-down 
from weak decays. 

Other models have also been used to study hadron ra- 
pidity distributions at RHIC. While results from the HI- 
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JING model |65|, Ill6l | are compatible with the observed 
charged particle rapidity density, the model does not in- 
clude interactions among minijet partons and final-state 
interactions among hadrons. The saturation model with- 
out final-state interactions also reproduces the experi- 
mental data [3^ Ib^ . Furthermore, the hadronic cas- 
cade model LUCIFER predicts a charged particle 
multiplicity at mid-rapidity that is comparable to the 
RHIC data On the other hand, the LEXUS model 
|l7|. which is based on a linear extrapolation of ultra- 
relativistic nucleon-nucleon scattering to nucle us- nu cleus 
collisions, predicts too many charged particles |ll7j com- 
pared with the PHOBOS data. The URQMD model [3 
also failed 114] in describing the charged particle multi- 
plicity at RHIC. 
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FIG. 21: (Color online) Rapidity distributions of charged par- 
ticles in central (& < 3 fm) Au-|-Au collisions at ^/sNN ~ 200 
GeV from the default AMPT model (solid curves), with- 
out final-state interactions (long-dashed curves), without nu- 
clear shadowing (dashed curves), and with jet quenching of 
dE/dx — —1 GeV/fm (dotted curves). 



To see the effect of hadronic interactions, we show in 
Fig. 1211 by long-dashed curves the rapidity distributions 
of charged particles obtained from the default AMPT 
model without final-state interactions (i.e., without par- 
ton and hadron cascades) for central (fe < 3 fm) Au+Au 
collisions at ^snn = 200 GeV. Compared to the case 
with hadronic scatterings, there is a significant increase 
in the numbers of total charged particles and pions at 
mid-rapidity. The kaon number, on the other hand, 
only increases slightly when production from final-state 
hadronic interactions is excluded. The ratios of p/p and 
K~^/tt~^ at mid-rapidity are 0.85 and 0.13, respectively, 
in the absence of final-state interactions, instead of 0.81 
and 0.17 from the default AMPT model including the 



hadron cascade. We note that although the default HI- 
JING "gF, 'IIG] with original a and b parameters gives a 
total charged particle multiplicity at mid-rapidity that is 
consistent with the RHIC data, including hadronic scat- 
terings reduces appreciably the final number. 

We also find that excluding parton cascade in the 
AMPT model changes the final charged particle yield at 
mid-rapidity at ^/sNN = 200 GeV by less than 5%. This 
indicates that hadron yields are not very sensitive to par- 
ton elastic scatterings in the default AMPT model. To 
take into account the effect of parton inelastic collisions 
such as gluon radiation, which is not included in ZPC 
as it includes at present only elastic scatterings, we have 
included in the AMPT model the default HIJING jet 
quenching, i.e., an energy loss of dE/dx = —1 GeV/fm, 
before minijet partons enter the ZPC parton cascade. Re- 
sults obtained with jet quenching for central Au-f Au col- 
lisions at y^SNN = 200 GeV are shown in Fig. |^ by 
dotted curves. We see that the quenching effects are 
larger for pions than for kaons, protons, and antiprotons. 
Since present calculations from the AMPT model with- 
out jet quenching already reproduce the data at collision 
energy of y^SNN = 200 GeV, and further inclusion of jet 
quenching of dE/dx = — 1 GeV/fm increases the final 
yield of total charged particles at mid-rapidity by about 
14%, our results for the rapidity distribution of charged 
particles are thus consistent with a weak jet quenching at 
this energy. However, the dense matter created in heavy 
ion collisions expands rapidly, thus the energy loss at the 
early stage may still be large |ll8l Ill9 ^ . 

The effect of initial nuclear shadowing of the parton 
distributions in nuclei is also shown in Fig. [2] Without 
initial nuclear shadowing effect on minijet parton produc- 
tion, the charged particle multiplicity at mid-rapidity in 
central Au-|-Au collisions at ^snn — 200 GeV increases 
by about 23%. Although this increase can be offset by 
using different parameters in the Lund string fragmenta- 
tion, to reproduce both SPS and RHIC data with same 
parameters in the default AMPT model requires the in- 
clusion of nuclear shadowing on minijet parton produc- 
tion. 



B. Particle ratios 

The energy dependence of the yields of particles and 
their ratios at mid-rapidity from central Au-I-Au colli- 
sions at AGS to central Au-f Au collisions at RHIC en- 
ergies are shown in Fig. [221 Curves represent AMPT re- 
sults, while data from AGS, SPS and RHIC experiments 
are shown by symbols: triangles for pions, squares for 
kaons and circles for protons and antiprotons in the up- 
per figure. The AGS data are from the E917 and E866 
Collaborations for central Au-f Au collisions at labora- 
tory kinetic energies of 7.94 and 10.7 GeV per nucleon 
|120( . from the E895 Collaboration for central Au-fAu 
collision s at the laboratory kinetic energy of 8 GeV per 
nucleon |l2lj . and from the E802 Collaboration for cen- 
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C. Baryon and antibaryon production 
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FIG. 22: (Color online) Energy dependence of particle yields 
and ratios of K~ /K^ , p/p, and /ir'^ at mid-rapidity in 
central Au+Au or Pb+Pb collisions. 



tral Au-KAu collisions at Piab = 11.6—11.7 GeV/c per 
nucleon |l22l |. The SPS data are from NA49 for central 
Pb-|-Pb collisions at laboratory energies of 40A, 80^ and 
158A GeV .105J. T he RHIC data are from the PHENIX 
Collaboration |l23l [T2i | and the STAR Collaboration 
|l25l Il26| for central Au-|-Au collisions at ^/snn — 130 
and 200 GeV. 

Since the antiproton yield increases with energy al- 
most logarithmically while the proton yield initially de- 
creases at low energies, the p/p ratio (data represented 
by circles) increases rapidly from almost at AGS to 
about 0.1 at the SPS energy of 158^ GeV then rapidly 
to about 0.8 at the maximum RHIC energy of ^/snn — 
200 GeV, indicating the formation of a nearly baryon- 
antibaryon symmetric matter at RHIC. The K'^ /n^ ra- 
tio (data represented by triangles), on the other hand, 
is almost constant from the SPS energy of 158A GeV 
to RHIC, suggesting an approximate chemical equilib- 
rium of str ange ness in heavy ion collisions in this en- 
ergy range |l27j . The /K^ ratio (data represented 
by squares) increases from below 0.2 at AGS to about 
0.6 at the SPS energy of 158^1 GeV, then gradually to 
about 0.9 at ^snn = 200 GeV, and the value near 1 at 
the top RHIC energy is closely related to the fact that 
the matter formed at RHIC is nearly baryon-antibaryon 
symmetric. 



FIG. 23: (Color online) Net-proton and antiproton rapidity 
distributions from the default AMPT model for central (6 < 3 
fm) Au-j-Au collisions at RHIC energies of y^sjvjv = 200 GeV 
(solid curves) and 56 GeV (dashed curves) versus data from 
the BRAHMS Collaboration for 5% most central collisions at 
= 200 GeV. 



Fig. 1231 shows the net-proton and antiproton rapidity 
distributions from the default AMPT model at two differ- 
ent RHIC energies. Filled and open circles are BRAHMS 
data on net-proton and antiprotons without taking into 
account corrections due to weak decay of hyperons 12^, 
i.e., the data include feed-down from hyperon weak de- 
cays. Since the AMPT results on proton and antiprotons 
exclude feed-down from weak decays, the antiproton yield 
at mid-rapidity is actually over-predicted by the AMPT 
model, as can be more clearly seen in Fig. 1201 while the 
net-proton rapidity distributions may agree reasonably 
with the BRAHMS data. 

It is further seen that the antiproton yield at mid- 
rapidity increases rapidly with collision energy, and peaks 
of net-proton and proton distributions shift toward larger 
rapidity at higher collision energies. Since the proton 
yield at mid-rapidity first decreases from the AGS en- 
ergy to the RHIC energy of ^Jsnn — 56 GeV and even- 
tually increases slowly with energy (see Fig. 122(1 . it may 
have a minimum between the SPS energy and the max- 
imum RHIC energy. On the other hand, the net-proton 
yield at mid-rapidity decreases with the colliding energy 
in the energy range from the AGS energy to the maxi- 
mum RHIC energy. 
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D. Transverse momentum spectra 




^ p, (GeV/c) 

FIG. 24: (Color online) Transverse momentum spectra of 
mid-rapidity pions, kaons, protons and antiprotons from cen- 
tral (6 < 3 fm) Au-|-Au collisions at y^sjvjv = 200 GeV from 
the default AMPT model. Data are from the PHENIX Col- 
laboration. 

For transverse momentum spectra, results from the de- 
fault AMPT model for pions, kaons, and protons at mid- 
rapidity from central (6 < 3 fm) Au+Au collisions at 
^/sWN = 200 GeV are shown in Fig.l24ltogether with the 
5% m ost central data from the PHENIX Collaboration 
[l2^ . Below pT = 2 GeV/c, the default AMPT model 
gives a reasonable description of the pion and kaon spec- 
tra It over-predicts, however, both the proton and 
antiproton spectra at lowpx- We note that the PHENIX 
data have been corrected for weak decays. 

The effect due to final-state hadronic scatterings on 
pion and proton transverse momentum spectra is illus- 
trated in Fig. 1251 where solid and dashed curves with 
statistical errors are, respectively, the results with and 
without hadron cascade in the default AMPT model. 
As found in Fig. ^| at SPS energies, hadronic rescat- 
terings increase significantly the inverse slope of the pro- 
ton transverse momentum spectrum while they do not 
affect much that of pions. As a result, the final pro- 
ton yield at mid-rapidity becomes close to the pion yield 
at jPj ~ 2 GeV/c, as observed in experiments at RHIC 
|l23l Il29j| . Without final-state scatterings, the proton 
yield given by the dashed curve is well below the pion 
yield up to ~ 1-7 GeV/c when statistical fluctua- 
tions in the AMPT calculations become large. Results 
from the default AMPT model thus indicates that the 
observed large p/n ratio at px 2 GcV/c is due to the 
collective transverse flow generated by flnal-state rescat- 




FIG. 25: (Color online) Transverse momentum spectra of 
mid-rapidity pions and protons from the default AMPT model 
with (solid curves) and without (dashed curves) final-state in- 
teractions in central Au-|-Au collisions at ^/snn = 200 GeV. 



terings. 

The above results are obtained without quenching of 
minijet partons due to gluon radiations. Including this 
effect would reduce the yield of hadrons at large trans- 
verse momentum as shown in Fig. 1261 leading to an ap- 
preciable discrepancy between the AMPT results and the 
experimental data from the PHENIX collaboration for 
hadrons with momenta greater than 1.5 GeV/c. This 
discrepancy may stem, however, from the fact that the 
AMPT model already underpredicts the STAR data for 
hadron yield above 1 GeV/c in pp collisions as shown in 
Fig. Since the STAR results are consistent with pre- 
vious UA5/CDF measurements at similar multiplicities 
'126^ but have signiflcantly higher mean kaon and pro- 
ton transverse momenta than interpolated values from 
the UA5/CDF pp/pp data shown in Fig. ^ the STAR 
NSD trigger might be quite different from the UA5/CDF 
NSD trigger. To make a more definitive conclusion on 
jet quenching in AA collisions thus requires a focused 
and higher-statistics study of pt spectra in both pp and 
Au+Au colhsions. 



E. Phi meson yield via dilepton and dikaon spectra 

Since the (f) meson is unstable, it can only be detected 
from its decay products such as the kaon-antikaon pair or 
the lepton pair. In the AMPT model, we follow the pro- 
duction and decay of (j) mesons as well as the scattering 
of the kaon daughters from their decays. Since the lat- 
ter destroys the possibility of reconstructing the parent 



18 




^ p, (GeV/c) 

FIG. 26: (Color online) Transverse momentum spectra with 
and without quenching of dE/dx — —1 GeV/fm for cen- 
tral Au+Au collisions at y^sjvjv ~ 200 GeV from the default 
AMPT model versus data from the PHENIX Collaboration. 




FIG. 27: (Color online) Transverse momentum spectra of 
mid-rapidity pions, kaons and protons for pp collisions at 200 
GeV from the default AMPT model versus data from the 
STAR Collaboration. 

(j) mesons, only cj) mesons with their decay kaons not un- 
dergoing scattering can be reconstructed from the kaon- 
antikaon invariant mass distribution. The kaon-antikaon 
pair from a ip meson decay is likely to undergo rescatter- 
ings in the medium, which would lead to a reconstructed 
invariant mass outside the original cj) meson peak. Hence 



(j) mesons decaying in the dense medium are difficult 
to be identified via reconstructed kaon-antikaon pairs. 
In contrast, dileptons have negligible final-state interac- 
tions with the surrounding hadronic medium, they are 
thus considered to carry useful information about hadron 
properties in hot dense hadronic matter |il3Q . il31j |. which 
are expec ted to be different from those in free space 

Since dimuons are emitted continuously during evolu- 
tion of the system, the total number of dimuons is given 

by 

Jo 

+ A^4W) ^^7^^" , (17) 
l0 

where Nip(t) denotes the number of (j) mesons at time t 
and T^^^+^-{M) is its decay width to dimuon as given 
by 

w.-(Ai)=c,.,-^(^i-^j [y+^ym 

The coefficient C^+^- = aV^^Cg^^i?/'^'^) = 1-634 x 
10~^ is determined from measured width of free (j) 
son at invariant mass M = m^. In Eq. (|17|l . the first 
term on the right-hand-side corresponds to dimuon pro- 
duction before a cut-off time tcut during the hadronic 
phase of heavy ion collisions, while the second term refers 
to dimuon emission after that cut-off time. The recon- 
structed (f) meson number is then obtained by dividing 
the above expression by the dimuon branching ratio in 
free-space of r^^fj^+^- /T^ — 3.7 x 10""*. The number of 
kaon-antikaon pairs coming from cj) meson decays can be 
similarly expressed as Eq. itTTI) : 

Nkk^ r^dt iV^(i)r^_^;^(M)-f 7V^(t,„t)%^-(19) 
Jo -L 

The (j) meson number from kaon-antikaon decays is then 
obtained by dividing Eq. H19() by the KK branching ratio 
in free-space r^^^^^/T^. 

Using the default AMPT model, we have studied the 
(j) meson yield reconstructed from K~^K~ and fJ-^fi~ 
pairs for central Pb-|-Pb collisions at SPS and for cen- 
tral Au-I-Au collisions at RHIC both cases we 
have found that, due to hadronic rescattering and ab- 
sorption in the kaonic channel, the (f> meson abundance 
at mid-rapidity is larger in the dimuon channel than in 
the dikaon channel, and the inverse slope parameter ob- 
tained from the transverse mass spectra of mesons in 
the range < tot — "i^ < 1 GeV/c^ is larger in the 
dimuon channel than in the dikaon channel. These fea- 
tures are consistent with the data at SPS from the NA49 
|134| and NA50 135] collaborations. Comparison of the 
results for RHIC with future experimental data will al- 
low us to learn if enhanced (j) meson production is due 
to the formation of the quark-gluon plasma during the 
early stage of collisions |47|. 
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the QGP is formed in heavy ion colhsions. In particu- 
lar, the observed abnormal suppression of J /if) in cen- 
tral Pb-|-Pb collisions at SPS has been attributed to the 
formation of the quark-gluon plasma in these collisions 
However, there are other possible mechanisms 
for J /ip suppression in heavy ion collisions 152L.153,. 1541]. 
e.g., J/?/; may be destroyed by collisions with incoming 
nucleons or with gluons in the initial partonic matter 
[155 , 1561 or with comoying hadr ons hi the hadr onic mat- 
ter liH^liElEsllTeSITMl^ 



F. Multistrange baryon production 

One possible signal for the quark-gluon plasma is en- 
hanced production of strange particles ^136j , particularly 
those consisting of multistrange quarks such as 2 and f2 
baryons as well as their antiparticles. The argument is 
that the rate for strange hadron production is small in 
hadronic matter due to large threshold and small cross 
sections while the production rate for strange quarks is 
large in the quark-gluon plasma ■ 

A detailed study of multi-strange baryons in the 
AMPT model Q shows that, although few multi- 
strange baryons are produced in HIJING (see Table I 
of Ref. ^3)' including the strangeness-exchange inter- 
actions hsted in Sec. HID 4l in the default AMPT model 
leads to an enhanced production of multistrange baryons 
in heavy ion collisions at both SPS and RHIC. For cen- 
tral Pb-|-Pb collisions at SPS, the strangeness-exchange 
reactions enhance modestly the yields of A, E, and S, 
but increase the O yield by more than an order of mag- 
nitude. However, the yield in central Pb-|-Pb collisions 
at SPS from the default AMPT model is still about a 
factor of 2 lower than the data, and this may indicate 
that strangeness production is already enhanced during 
early stage of collisions. 

Predictions from the default AMPT model Q also 
show that the slope parameters obtained from the trans- 
verse mass spectra of multistrange baryons reveal a 
plateau structure since these particles, mostly gener- 
ated by strangeness-exchange reactions in the model, are 
weakly interacting and decouples rather early from the 
system. It will be very interesting to test these predic- 
tions against the RHI C data. In particular, the quark 
coalescence model [lH [Ha 140. ML ImIImI 
Il45| . which assumes that all energy in the early stage of 
RHIC collisions is in the partonic degrees of freedom as 
assumed in the string melting scenario 4^ n^Sn of the 
AMPT model, relates the transverse momentum spectra 
or the elliptic flow of multistrange baryons with other 
hadrons such as kaons and protons. These relations can 
be drastically different from the default AMPT model, 
where the partonic stage includes only minijet partons 
and has a much smaller effect than later hadronic stage. 

G. Equilibration between J/iIj and open charm 



llMll69|,lMlIl3- 

For heavy ion collisions at RHIC energies, multiple 
pairs of charm- ant icharm quarks can be produ ced i n one 
event. Estimates using the kinetic approach (17^ have 
shown that the number of J /tp produced from the inter- 
action between charm and anticharm quarks, i.e., the in- 
verse reaction of J /ip dissociation by gluons, may exceed 
that expected from the superposition of initial nucleon- 
nucleon i nteractions. C alculations based on the statisti- 
cal model |l73lll74lll75l |. which assumes that J/ip formed 
during hadronization of the QGP is in chemical equilib- 
rium with charm mesons, also predict that the J num- 
ber is comparable to the expected primary yield. 

To study the above discussed effects on J/ij] produc- 
tion, the default AMPT model has been modified |176| 
to include J /ip absorption and production in both ini- 
tial partonic and flnal hadronic matters |4ll |. It was 
found that because of these final-state interactions the 
final J/ ijj yield in central Au+ Au collisions at RHIC may 
exceed the suppressed primary J yield which is based 
only on the color screening mechanism. Furthermore, the 
J/f/; yield at RHIC depends on the charm quark mass in 
the partonic matter and charmed hadron masses in the 
hadronic matter. 

We note that the elliptic flows of J /'ip and charmed 
hadrons are related in the quark coalescence model, and 
they exhibit novel mass effects when the constituent 
quark masses in the hadron are different |l42j | . The trans- 
verse momentum spectra of electrons from open charm 
deca ys ar e, however, insensitive to the charm transverse 
flow |l77j |. Com paring RHIC data on charm hadron or 
electron V2 jl43i | with results from the quark coalescence 
model may allow us to determine whether in heavy ion 
collisions charm quarks flow collectively as light quarks. 



J/tp production has long been proposed as a possible 
signal for the formation of t he qu ark-gluon plasma in rel- 
ativistic heavy ion collisions 146]. Based on the expecta- 
tion that the effective potential between charm and an- 
ticharm quarks changes in the QGP due to color screen- 
ing effect, they will not form bound states above certain 
critical temperature, which is somewhat higher than the 
deconfinement temperature. The qualitative change in 
the heavy quark effective potential h as been verifie d by 
results from lattice QCD simulations |l47lll48l ITigl] . As 
a result, J/ip production is expected to be suppressed if 



H. Elliptic flow 

Elliptic fiow in heavy ion collisions is a measure of 
the asymmetry of particle momentum distributions in 
the transverse plane and is generated by the anisotropic 
pressure gradient in initial hot dense matter as a re- 
sult of the spatial asyrnmetr y in non-central collisions 
[m [Izl [Ha nil [isl [iS^. It is defined as one half 
of the second Fourier coefficient of the azimuthal angle 
distribution of particle transverse momentum and can be 
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evaluated as 



V2 




(20) 



In the above, the x-axis is along the impact parameter in 
the transverse plane of each event while the y-axis points 
out of the reaction plane, and the average is taken over 
particles in consideration. 

At RHIC, large hadron elliptic flows have been ob- 
served Hm, 185, ,186, 187, 188, 189, IgOL and their 
dependence on transvers e momentum jl84l Il85l | and 
pseudorapidity |l86l Il9lj as well as on particle species 
[l87, 192] have also been studied. To understand these 
experimental results, many theoretical models have been 
introduced, and these include semi-an alytic mod els |l93l | , 
models based on pa rton^ e nergy loss |194 Il95j , hydro- 
dynamic models 9, 180, 196, '191], transport models 
El mi [m JM, m ml, including a hybrid model 
that uses the transpo rt m odel as an after-burner of the 
hydrodynamic model '202'], and quark coalescence mod- 
els [m^ J38, 139, 140, 141, 142, 143, 144, 145j. Trans- 
port models based on hadronic and/or string degr ees of 
freedom in general give a smaller elliptic flow |187| than 
that observed at RHIC. Although hydrodynamic mod- 
els can explain the large elliptic flow at low transverse 
momenta with a large initial energy density, they pre- 
dict a continuously increasing elliptic flow with increas- 
ing hadron transverse momenta instead of the observed 
level off in experiments 187, 189. 203. 204j , indicating 
that high transverse momentum particles do not reach 
thermal equilibrium. With the AMPT model, we can 
address such non-equilibrium effects to obtain informa- 
tion on the degree of thermalization in these collisions. It 
further allows us to study the effects due to both partonic 
and hadronic rescatterings. 

In Fig.|2Hl we show the elliptic flow of all charged parti- 
cles in Au-|-Au collisions at ^/sNN = 200 GeV from both 
default AMPT model and the extended model with string 
melting as a function of A'part , which is the total number 
of participant nucleons after primary collisions but be- 
fore partonic and hadronic rescatterings |49j| . To compare 
with PHOBOS data l205l|, we have included only charged 
particles with pseudorapidities and transverse momenta 
in the ranges 77 G (—1,1) and px & (0.1,4) GeV/c, re- 
spectively, for evaluating the elliptic flow. Error bars 
in the figure represent the statistical error in our calcu- 
lations. We see that the value of V2 depends not only 
on whether initial strings are converted to partons but 
also on the parton cross section used in the model, sim- 
ilar to that seen in heavy ion collisions at ^/snn = 130 
GeV 45]. With the default AMPT model, which includes 
only scattering among minijet partons, the elliptic flow 
obtained with a parton scattering cross section of 3 mb 
is found to be too small to account for the experimental 
data from the PHOBOS collaboration. The elliptic flow 
is larger in the string melting scenario as soft partons 
from melted strings also participate in partonic scatter- 
ings, and its value is further increased when the parton 
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FIG. 28: (Color online) Centrality dependence of charged 
hadron elliptic flow in Au-|-Au collisions at ^/s nn = 200 
GeV. Data from the PHOBOS collaboration 205] are shown 
by circles, while results from the AMPT mode with diflPerent 
partonic dynamics are shown by curves. 



scattering cross section increases from 3 mb to 10 mb, 
indicating that the dense system created in heavy ion 
collisions at RHIC does not quite reach thermal equilib- 
rium. 

In Figure EHl the differential elliptic flows V2{pt) of 
pions, kaons and protons from the AMPT model in the 
string melting scenario with a parton scattering cross sec- 
tion of 6 mb are shown for minimum-bias Au+Au colli- 
sions at y/sWw = 200 GeV. Although the PHENIX data 
are for particles with pseudorapidities i] G (—0.35,0.35), 
the AMPT results are for 77 g (—1, 1) in order to improve 
the statistics of our calculation. In this low pt region, we 
observe the mass ordering of W2(pt), i-e., hadrons with 
larger masses have smaller V2 values at the same pj , sim- 
ilar to that from the hydrodynamic model 0, Q. For 
hadrons at higher px, which are not studied here, a scal- 
ing of hadron elliptic flows according to their constituent 
quark content ha s been proposed bas ed on the quark co- 
alescence model |l37l Il38l Il39l |206j. Except for pions, 
this scaling seem s to b e confirmed by experimental data 
from both STAR |203 and PHENIX ^ 2o|] coUaborations. 
The violation of pion elliptic flow from the constituent 
quark number scaling ha s bee n attributed to effects due 
to reson ance decays |209l l2icij and the binding energy of 
hadrons Il42ll. 



The pseudorapidity dependence of the elliptic flow in 
minimum-bias Au-j-Au collisions at y^s^N = 2 00 GeV 
has also been studied in the AMPT model |21lj . While 
the string melting scenario with a parton scattering cross 
section of 3 mb as well as the default AMPT model gives 
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FIG. 29: (Color online) Transverse momentum dependence of 
elliptic flow of identified hadrons for minimum-bias Au+Au 
collisions at ^/sNN ~ 200 GeV. Circles are data from the 
PHENIX Collaboration .2Q8,] while curves are results from 
the AMPT model with string melting and using 6 mb for the 
parton scattering cross section. 



a better description of the PHOBOS data | 205| at large 
pseudorapidity, the string melting scenario with a larger 
parton scattering cross section reproduces the PHOBOS 
data at mid-pseudorapidity. This may not be unreason- 
able as not all strings are expected to melt at large pseu- 
dorapidity where the smaller particle multiplicity leads 
to a lower energy density. 

The AMPT model has further been used to study 
higher-order anisotropic flows and vq of charged 
hadrons at mid-rapidity in heavy ion collisions at RHIC 
||212J. It was found that the same large parton scat- 
tering cross section used in explaining the measured V2 
of charged hadrons could also reproduce recent data on 
their V4 an d from Au-|-Au collisions at ^snn — 200 
GeV |21.l| |. Furthermore, the V4 was seen to be a 
more sensitive probe of the initial partonic dynamics in 
these collisions than V2- Moreover, higher-order parton 
anisotropic flows are non-negligible and satisfy the scal- 
ing relation Vn.q{pT) ^ '^2^q(j>T)j which leads naturally 
to the observed similar scaling relation among hadron 
anisotropic flows when the coalescence model is used to 
describe hadron production from the partonic matter. 



not only on the spatial exten t of the emission so urce but 
also on its emission duration |215l l216l l217l |218| | . In par- 
ticular, the long emission time as a result of the phase 
transition from the quark-gluon plasma to the hadronic 
matter in relativistic heavy ion collisions may lead to 
an emission sou rce which has a i?out/^sidc ratio much 
larger than one |218l l219l l22Clj , where the out-direction 
is along the total transverse momentum of detected two 
particles and the side-direction is perpendicular to both 
the out-dir ection an d the beam direction (called the long- 
direction) mini. Since the quark-gluon plasma is ex- 
pected to be formed in heavy ion collisions at RHIC, it is 
thus surprising to find that the extracted ratio i?out/^^sidc 
from a Gaussian fit to the measured correlation function 
of two identic al pions in Au+Au collisions is close to one 
■22?, "22?, "2235. This is in sharp contrast with calculations 
from hydrodynamic models |219l |224) , where i?out /-Rsidc 
is typically well above one. 

Denoting the single-particle emission function by 
5(0;, p), the HBT correlation function for two identical 
bos ons i n th e absence of final-state interactions is given 
by IMIll: 



C2(Q,K) = 1 + 

/ d'^x\d'^X2S{xi,J^)S{x2,^) cos [Q-[xi — X2)\ 

J d'^XiS{xi,Pi) J d'^X2S{x2,P2) 



,(21) 



where pi and P2 are momenta of the two hadrons, K — 
(Pi + P2)/2, and g = (pi - p2,£;i - E2). The three- 
dimensional correlation function in Q can be shown as 
one-dimensional functions of the p rojections of Q in the 
"out-side- long" coordinate system |216l l217j . 

Expecting that the emission function is sufhciently 
smooth in momentum space, the size of the emission 
source can then be related to the emission function as: 



ia2C2(Q,K) 



2 dQ.dQ, 



Q=0 



These source radii are thus expressed in terms of the 
space-time variances, D^^y = {x ■ y) — {x){y), with (x) 
denoting the average value of x. In the above, Xi's {i = 
1 — 3) denote the projections of the particle position at 
freeze-out in the "out-side- long" system, i.e., a;out, Xsidc 
and xiong, respectively; and (3 = J^/Kq with Kq being 
the average energy of the two particles. 

The experimentally measured two-particle correlation 
function C2(Q, K) in central heavy ion collisions is usu- 
ally fitted by a four-parameter Gaussian function after 
correcting for final-state Coulomb interactions, i.e., 



I. Two-particle interferometry 



C2(Q, K) = 1 + A cxp - ^ Rl {K)Ql 



(23) 



The Hanbury-Brown Twiss (HBT) effect was first used 
to measure the size of an emission source like a star |214j . 
For heavy ion collisions, HBT may provide information 



If the emission source is Gaussian in space-time, the fit- 
ted radii Ra would be identical to the source radii deter- 
mined from the emission function via Eq. (I22II . However, 
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because of space-time correlations in the emission func- 
tion, such as those induced by the collective flow, the 
fitted radii can be quite different from the source radii 
[22^11227^ . 

We have used the AMPT model to study the in- 
terferometry of two identical pions or kaons in central 
Au+ Au collisions at RHIC 0, |H3| • The source of the 
emitted particles have been obtained from their space- 
time coordinates x and momenta p at kinetic freeze-out, 
i.e., at their last interactions. The correlation function 
C2(Q, K) is then evaluated in the frame of longitudinally 
comoving system using the program Correlation After 
Burner |228| . In these calculations, the cut-off time tcut 
for hadron cascade has been chosen as 200 fm/c for HBT 
studies instead of the default 30 fm/c as we are inter- 
ested in the space-time and momentum distributions of 
hadrons at freeze-out even though their rapidity distri- 
butions and momentum spectra essentially do not change 
after 30 fm/c. 
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FIG. 30: (Color online) Pion source radii with (dashed curves) 
and without (solid curves) uu decays as well as the fitted radius 
and A parameters from the Gaussian fit to the correlation 
function (curves with squares) at mid-rapidity as functions 
of fcx from AMPT with string melting and ap — 6mb for 
central Au+Au collisions at ^ysNN = 200 GeV. Circles are 
the PHOBOS data, while diamonds are the STAR data. 

For central {b = fm) Au-|-Au collisions at ^/snn = 
200 GeV, results for charged pions from the AMPT model 
with string melting and parton cross section dp = 6 mb 
are shown in Fig. 1301 together with experimental data 
from the STAR ^] and PHOBOS 230J Collaborations. 
The STAR data at the three transverse momentum /ct 
bins of (0.15,0.25), (0.25,0.35), and (0.35,0.60) GeV/c 
are for mid-rapidity and 0-5% most central collisions, 
while the PHOBOS data are for 0-15% most central col- 
lisions with {y-wTr) = 0.9. We find that the source radii 
-Rout (upper-left panel), i?sido (upper-right panel), and 
^long (lower-left panel) of pions including those from 



uj decays (dashed curves) are a factor of 2 to 3 larger 
than the radius parameters from a Gaussian fit to the 
three-dimensional correlation function obtained from the 
AMPT model (curves with squares) or to the experimen- 
tally measured one (open or filled circles). Excluding pi- 
ons from Lu decays reduces the source radii (solid curves) , 
and this brings i?sidc close to the fitted one while i?out 
and i?iong can still be a factor of 2 larger than fitted 
ones. The emission source from the AMPT model thus 
deviates appreciably from a Gaussian one, as found previ- 
ously instudies at SPS using the RQMD transport model 
226. l 22'i| . In this case, it will be useful to compare the 
emission source from the AMPT model with that ex- 
tracted from measured two -part icle correlation functions 
using the imaging method |23lL l2^ 
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FIG. 31: (Color online) Radii, A parameters and -Rout/^^sido 
for mid-rapidity pions with 125 < px < 225 MeV/c from 
Gaussian fits of the correlation function as functions of ap. 
The STAR data at ^/sj^ = 130 GeV and 200 GeV are shown 
by open and filled symbols, respectively. 



Since Eq. gives 



and i?sid 



D 



- 2 + Dp^t,p^u (24) 

the ratio -Rout/^sido contains 



information about the du ration of emiss i on a nd has thus 
been studied extensively Ull IHlllll IllH . However, 
the Xout — t distributions at freeze-out from the AMPT 
model show a strong positive Xout — t correlation for 
both pions and kaons [IsL Isoj . This leads to a positive 
■DxoMt,P±t and thus a negative contribution to R^^^ that 
can be as large as the positive duration-time term, mak- 
ing it difficult to extract information about the duration 
of emission from the ratio i?out /^f^side alone. 

The dependence of fitted radii on cTp, with CTj, = de- 
noting the default AMPT model without string melting. 
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is shown in Fig. 1311 for mid-rapidity charged pions with 
125 < pt < 225 MeV/c. Also shown by circles, trian- 
gles, diamonds, and squares are the STAR data from 
central Au-fAu collisions for mid-rapidity pions with 
125 < PT < 225 MeV/c at V^jvw = 130 GeV (open 
symbols) and for pions with 150 < pt < 250 MeV/c at 
^/sWn = 200 GeV (filled symbols) for Rout, i?sido, ^long, 
and A parameter, respectively. Both the radius and the 
A parameters from the AMPT model with string melting 
and (Tp = 3—10 mb are close to the experimental values 



[22 1| . while results from the default AMPT model over- 
predict the A parameter. For results with string melting, 
all three radius parameters are seen to increase with in- 
creasing parton cross section ap presumably as a result 
of the larger source size at freeze-out, while the extracted 
A parameter decreases gradually with increasing <7p. 
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only dynamical model that gives a i?out /^sido ratio close 
to 1 and also gives roughly the correct magnitude for the 
fitted radii. Future studies using this transport model 
on the X — t correlation term |233l | or on azimuthal HBT 
in non-central heavy ion collisions will help us to further 
test the AMPT model and understand the freeze-out dy- 



VI. RESULTS AT LHC 
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FIG. 33: (Color online) Rapidity distributions for central (6 < 
3 fm) Pb-|-Pb collisions at y^sjvjv ~ 5500 GeV from the default 
AMPT model with (curves with circles) and without (curves 
with no symbols) nuclear shadowing. 



FIG. 32: (Color online) Ront/Rside ratios for mid-rapidity pi- 
ons at ^/sNN ~ 200 GeV as functions of /ct . Solid and dashed 
curves are AMPT results from the emission function without 
and with pions from lu decays, respectively. The curve with 
stars is the AMPT result from the Gaussian fit to the two-pion 
correlation function. 

Fig. 1321 shows the ratio i?out/-Rsido for mid-rapidity 
charged pions from AMPT with string melting and 
(Tp = 6 mb for central (6 = fm) Au-|-Au collisions at 
y^sjvjv = 200 GeV. It is seen that the ratios obtained 
from the emission function, with solid and dashed curves 
for results without and with pions from lo decays, have 
values between 1.0 and 1.7, consistent with predictions 
from the hydrodynamic mode l wi th freeze-out treated 
via hadronic transport model |219| |. However, the ratio 
from the Gaussian fit to the three-dimensional correla- 
tion function (the curve with stars) is much closer to 1, 
similar to that extracted from a Gaussian fit to the mea- 
sured correlation function 221, 229, 230). 

The AMPT model with string melting is so far the 



The AMPT model can also be used to study heavy 
ion collisions at LHC. Fig. gives the charged parti- 
cle pseudorapidity distribution and the rapidity distribu- 
tions of charged pions, kaons, protons, and antiprotons 
with and without nuclear shadowing in central (& < 3 fm) 
Pb+Pb collisions at ^/snn — 5.5 TeV from the default 
AMPT model. The distributions are significantly wider 
and higher than corresponding distributions at RHIC. 
At mid-rapidity, the distributions without shadowing are 
higher than corresponding ones with shadowing by about 
80%. The highest value at mid-rapidity is about 4500, 
well within the LHC detector limit of 7000 particles per 
unit rapidity. The mid-rapidity density with nuclear 
shadowing is about 2500, more than a factor of three 
higher than that at RHIC. It is higher than the loga- 
rithmic extrapolation from lower energy data but lower 
than the saturation model prediction of about 3500 23'^ . 
The charged hadron pseudorapidity distribution shows a 
clear plateau structure which is very differ ent from pre- 
dictions from saturation models |235l l236l | . The proton 
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and antiproton rapidity distributions are close to each 
other and arc almost flat. This is different from the pro- 
ton and antiproton distributions at RHIC where protons 
clearly dominate at large rapidities. Note that the cut- 
off time tcut for hadron cascade has been chosen as 200 
fm/c at LHC instead of the default 30 fm/c due to the 
longer lifetime of the formed matter and the larger ra- 
pidity width in heavy ion collisions at LHC. 
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FIG. 34; (Color online) Transverse momentum spectra of pi- 
ons, kaons, and protons from the default AMPT model for 
central (6 < 3 fm) Pb-|-Pb collisions at ^snn = 5500 GeV 
(curves with circles) and for central (6 < 3 fm) Au-|-Au colli- 
sions at ^ysNN ~ 200 GeV (curves with no symbols). 



Transverse momentum spectra at LHC are shown in 
Fig-IMl It is seen that the inverse slope parameters, par- 
ticularly for kaons and protons with transverse momenta 
below 0.5 GeV/c and 1 GeV/c, respectively, are larger 
than at RHIC as a result of stronger transverse flows. 
Similar to that observed at RHIC, the proton spectrum 
is below that of pions at low transverse momenta, but 
they become comparable at about 2 GeV/c. As in heavy 
ion collisions at SPS (Fig. Uni) and RHIC (Fig.|23l, the 
strong transverse flow is due to final-state interactions, 
and this is shown in Fig. 1351 where the pion, kaon, and 
proton spectra at LHC obtained from the AMPT model 
with and without final-state interactions are compared. 

Because of the uncertainty in the parton distribution 
functions at small-x in heavy nuclei, the above results on 
LHC have large uncertainties. Currently, the initial con- 
dition of the AMPT model is obtained from the HIJING 
model with minijets and strings. The parton distribution 
functions used in the HIJING model are the Duke-Owens 
set 1 p2|, which are quite old. For example, Fig.l36lshows 
the gluon and u-quark distribution functions at — 9 
GeV^ from the Duke-Owens set 1 (circles), MRSA ^3 
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FIG. 35: (Color online) Transverse momentum spectra of pi- 
ons, kaons, and protons for central (6 < 3 fm) Pb-|-Pb colli- 
sions at y/s^ = 5500 GeV from the default AMPT model 
with (curves with circles) and without (curves with no sym- 
bols) final-state interactions. 



(squares), and the recent 4-flavor CTEQ6M |238l | (di- 
amonds) parameterizations. It is seen that the Duke- 
Owcns set 1 parameterization has far fewer partons at 
small-x, e.g., when x < 0.01, than the other two more 
recent parameterizations. At the top RHIC energy, a 
pair of 2 GeV back-to-back minijets at mid-rapidity cor- 
responds to X ~ 0.02 for the initial partons, and at the 
LHC energy it corresponds to cc ~ 0.00073. Thus small- 
X partons play much more important roles at LHC than 
at RHIC, and the AMPT model will have a much larger 
uncertainty at LHC due to the underestimate of small- a; 
partons from the Duke-Owens set 1. An update of the 
parton distribution functions for the HIJING model has 
been done 63] where the Gluck-Reya-Vogt parton distri- 
bution functions 239] have been implemented together 
with a new nuclear shadowing parametrization. To make 
predictions with better accuracy for LHC, we will need 
to update the parton distribution functions in a future 
version of the AMPT model. 



VII. DISCUSSIONS 

Two versions of the AMPT model have been used 
in the present study. The default AMPT model (cur- 
rent version 1.11), which includes only minijet partons in 
the parton cascade and uses the Lund string model for 
hadronization, is found to give a reasonable description of 
hadron rapidity distributions and transverse momentum 
spectra observed in heavy ion collisions at both SPS and 
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FIG. 36: (Color online) Gluon (curves with filled symbols) 
and u-quark (curves with open symbols) distribution func- 
tions (multiplied by x) at = 9 GeV^ from three different 
parameterizations. 
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FIG. 37: (Color online) Comparison of transverse momentum 
spectra of mid-rapidity pions from the AMPT model and the 
PHENIX data 123] for 5% most central Au-f An coUisions at 
Vsiviv = 200 GeV. 



RHIC. It under-predicts, however, the magnitude of the 
elliptic flow and also fails to reproduce the A parameter 
of the two-pion correlation function measured at RHIC. 
The latter can be described, on the other hand, by the 
AMPT model with string melting (current version 2.11) 
when the parton scattering cross section is about 6 mb as 
shown in Sects. IVTIl and lVTI This extended model under- 
predicts, however, the inverse slopes of hadron transverse 
momentum spectra, and also fails to describe the baryon 
rapidity distributions as shown in Fig. 1371 and Fig. 1381 
respectively. Since the AMPT model has not been able 
to describe all experimental observables at RHIC within 
a single version, further improvements are thus needed. 

The initial condition of the AMPT model is obtained 
from the HIJING model with minijets and strings. Al- 
though the string melting mechanism 45, 48, 50] is in- 
troduced to convert the energy in initial excited strings 
into partons in order to better model the partonic ini- 
tial condition at high energy densities, it is modeled by 
using the Lund string fragmentation to hadrons as an 
intermediate process and then converting these hadrons 
into their valence quarks. This is equivalent to staying 
at the parton level until strings in the Lund model have 
generated all the quark-antiquark pairs before forming 
hadrons. The current string melting can be viewed as 
a minimal implementation because, in the limit of no 
partonic interactions (and before hadron cascade starts), 
it reduces to HIJING results for all hadrons other than 
the flavor-diagonal mesons within the SU(2) flavor space, 
n^, 77, and lu. As a result, the initial partonic mat- 
ter in this scenario does not contain gluons. Although 



this is unphysical as one expects gluons to dominate the 
initial stage of ultra-relativistic heavy ion collisions, our 
study of the elliptic flow ^-^i^ the pion interferome- 
try |48l ] depends more on the effect of partonic scatter- 
ings instead of the composition of the partonic matter. 
However, for dilepton production from the partonic mat- 
ter, the flavo r composition of the partonic matter is im- 
portant 124(1 ]. and at present they cannot be addressed 
within the AMPT model. To extend the string melting 
scenario to include gluons in the AMPT model, we need 
to study the problem of both quark-antiquark and gluon 
production from a strong color field. Also, alternative 
descriptions of initial conditions, such as that from the 
parton saturation model [MIMIIlllllllllllllIll, 
may also be used as input to the AMPT model. 

The parton cascade model ZPC in the AMPT model 
only includes leading-order two-body partonic interac- 
tions (2-to-2), while higher order processes (m-to-n in 
general) , which might become dominant at high densities 
during the early stage of relativistic heavy ion collisions 
j24lj . have not been included. Unfortunately, it is dif- 
ficult to implement many-body interactions in a trans- 
port model, even for the next-leading-order processes 
{99 ^ 999) in * parton cascade [242]. Also, we currently 
treat the screening mass /i as a parameter for chang- 
ing the parton scattering cross sections. In principle, it 
should be evaluated dynamically from the effective tem- 
perature of the evolving parton system, leading thus to 
medium-dependent parton scattering cross sections. 

Furthermore, treating interactions at a distance as in 
the cascade model used for solving the Boltzmann cqua- 
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FIG. 38: (Color online) Rapidity distributions of total 
charged particles (upper left panel), protons and antiprotons 
(upper right panel), charged pions (lower left panel), and 
charged kaons (lower right panel) in 5% most central Au+Au 
collisions at ^/snn ~ 200 GeV from the default AMPT model 
(solid curves) and the AMPT model with string melting and 
CTp = 6 mb (dashed curves) . 



tions leads to possible violation of causality when the 
mean free path of partons is shorter than their interac- 
tion length, i.e.. 



A 



1 



< J- or > 1. 

pa \ TT n 



(25) 



If we take the parton density as p — 10/fm^ at RHIC and 
take the parton scattering cross section as tr = 3 mb, 
then we have p'^a^/n = 0.86 and thus expect that the 
usual method used in the cascade model for solving the 
Boltzmann equation may not be very accurate. However, 
the Boltzmann equation with only two-body scatterings 
is invariant under the transformation 



f ^ f I and a 



a 
T' 
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but the value of p^cr^ /tt is reduced by a factor of /. 
With sufficiently large I, the so-called parton subdivision 
can then overco me the causality pro blem in the cascade 
method HI El |20ll HI Hi |24|. Although parton 
subdivision has been implemented in ZPC, it remains to 
be implemented in the AMPT model mainly due to the 
complication in the hadronization process which converts 
partons into hadrons following the parton cascade. How 
the causality violation in the AMPT mode l affects the 
final observables in heavy ion collisions |20lj is yet to be 
studied. 

In the hadronization process in the default AMPT 
model, a minijet recombines with its parent string even 



though the minijet has gone through partonic interac- 
tions and thus has changed its color charge. The HIJING 
model also takes a similar approach in the jet quenching 
process because a gluon is still being associated with the 
same parent string after it goes through inelastic colli- 
sions. Our current prescription guarantees that fragmen- 
tation can successfully proceed for all strings after parton 
interactions, and it also enables the default AMPT model 
to reduce to HIJING results in the limit of no partonic 
(and hadronic) interactions. Although the freeze-out po- 
sitions and times of minijet partons are averaged in a 
string before its fragmentation, the current treatment 
of color configurations certainly needs to be improved, 
where new color-singlet strings over the whole volume 
need to be reconstructed after partonic interactions. 

To describe charged particle multiplicities in heavy ion 
collisions at SPS, we have changed the values of the a and 
h parameters in the Lund symmetric fragmentation func- 
tion given by Eq. Q based on the possibility that these 
parameters could be modified in the dense matter formed 
in heavy ion collisions. However, the AMPT model needs 
the default a and h values, i.e., the values in the HIJING 
model, in order to reproduce the charged particle multi- 
plicities in pp and pp collisions. The a and h values in the 
AMPT model is thus expected to depend on the atomic 
weight of colliding nuclei and the centrality of their col- 
lisions. Although the AMPT results on o? -f Au collisions 
using default a and h v alues k9|| agre e reasonably with 
the RHIC d + Au data pil |243, |24i|, the detailed de- 
pendence of a and b parameters on the system size has 
not been studied. 

At present, phase transition in the AMPT model with 
string melting is modeled by a simple quark coalescence 
using the current quark masses. The failure of the pro- 
ton rapidity distribution in this model, shown in Fig. 1381 
is related to this assumption as the proton mass is not 
generated dynamically but is given to the coalescing sys- 
tem of three light quarks. In the case that the invariant 
mass of the three quarks is small, the resulting proton 
with its physical mass tends to overpopulate at mid- 
rapidity. This problem may be avoided if protons are 
formed from coalescence of quarks with constituent quark 
masses. However, the current string melting mechanism 
would fail to convert strings into partons since the pion, 
as a Goldstone boson, cannot be decomposed into a quark 
and an antiquark with constituent quark masses. A more 
consistent method is perhaps to use density or temper- 
ature dependent quark masses, so that they correspond 
to current masses at high temperature and constituent 
masses at low temperature near the phase transition, 
where a scalar field is responsible for the changing masses. 
This method also has the advantage that it can qualita- 
tively describe the equation of state of the QGP jl82j |. 
Also, the current coalescence model could have prob- 
lems with entropy, because quark coalescence reduces the 
number of particles by a factor of two to three, although 
entropy also depends on the degeneracy and mass of pro- 
duced hadrons. We have studied the entropy problem 
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in a schematic model using a ther maliz ed QGP with the 
same parameters as given in Ref. |209| . Because of pro- 
ductions of massive resonances such as p, w, and A, the 
total entropy is reduced by about 15% when partons are 
converted to hadrons via coalescence, and this is mainly 
due to a similar violation in energy. To eliminate entropy 
violation thus requires a proper treatment of energy con- 
servation. This requires the inclusion of field energy in 
both QGP and hadronic matter, which is also necessary 
for properly describing the equation of state. Further- 
more, although the spatial correlation in the quark coa- 
lescence based on the closest neighbors leads to momen- 
tum correlation among coalescing partons in the pres- 
ence of collective flow, explicit momentum correlation 
is not considered in our quark coalescence model, and 
random momentum distributions of coalescing quarks 
tend to give the artificial peak at mid-rapidity shown 
in Fig. I2H1 The quark coalescence model can be im- 
proved by following the method used recently in studying 
hadron production from the quark-gluon plasma, which 
has been shown to describe satisfactorily the observed 
quark number scaling of hadron elliptic flows and large 
baryon to pion ratio at intermediate transv erse momenta 
[iHlISImlliilMlllilllilll^ Also, statis- 
tical models, which generate hadronic states according to 
their statistical weights, provides an alternative method 
to describe the parton-to-hadron phase transition. 

Many of the hadronic cross sections used in the 
hadronic stage of the AMPT model have not been studied 
theoretically in detail or well constrained by experimental 
data, and they may be important for some observables. 
For example, from the study at SPS energies with the 
default AMPT model, we have found that the explicit 
inclusion of K* mesons and baryon-antibaryon produc- 
tion from two-meson states are important for strangeness 
and antiproton production, respectively. In the AMPT 
model with string melting, the phase transition happens 
later than in the default AMPT model due to the larger 
energy and lifetime of the partonic system, hadronic ef- 
fects are less important than in the default AMPT model. 
Even in this case, it is essential to know the cross sections 
of rj meson interactions with other hadrons as they may 
affect the final 77 yield and consequently the height (or 
the A parameter) of the two-pion correlation function. 
Effects due to the uncertainties of these hadronic input 
parameters on final observables at RHIC will be studied 
in the future. 



VIII. SUMMARY 

To study high energy heavy ion collisions at SPS, 
RHIC and even higher energies such as at the LHC, a 
multi-phase transport (AMPT) model has been devel- 
oped. It consists mainly of four components: the HI- 
JING model to convert the initial incident energy to 
the production of hard minijet partons and soft strings, 
with excited strings further converting to partons in the 



AMPT model with string melting; Zhang's parton cas- 
cade (ZPC) to model the interactions among partons; 
the Lund string fragmentation as implemented in JET- 
SET/PYTHIA to convert the excited strings to hadrons 
in the default model or a simple quark coalescence model 
to convert partons into hadrons in the case of string melt- 
ing; and the extended relativistic transport (ART) model 
for describing interactions among hadrons. In this pa- 
per, we have described in detail the physics input in each 
component and how the different components are com- 
bined to give a comprehensive description of the dynam- 
ics of relativistic heavy ion collisions. We have used this 
model to study various observables in heavy ion colli- 
sions and to address the relative importance of partonic 
and hadronic effects on these observables. In particu- 
lar, the AMPT model has been used to study the rapid- 
ity distributions of particles such as pions, kaons, pro- 
tons and antiprotons, their transverse momentum spec- 
tra, the elliptic flow, and the interferometry of two iden- 
tical mesons. We find that the default AMPT model 
(current version 1.11) gives a reasonable description of 
rapidity distributions and transverse momentum spectra, 
while the AMPT model with string melting (current ver- 
sion 2.11) describes both the magnitude of the elliptic 
flow at mid-rapidity and the pion correlation function 
with a parton cross section of about 6 mb. 



High energy heavy ion collisions is a complex process 
involving the initial conditions, the interactions of ini- 
tially produced partons and of later hadronic matter as 
well as the transition between these two phases of matter. 
The AMPT model is an attempt to incorporate these dif- 
ferent physics as much as we can at present. Since there 
are many uncertainties in the AMPT model, the model is 
currently more a simulation tool rather than a finalized 
code. Nevertheless, this model provides the possibility 
to study the dependence of various observables on these 
physical effects. For example, we have found that the 
elliptic flows are sensitive to early parton dynamics, and 
the HBT interferometry is instead affected by the com- 
plicated late hadron freeze-out dynamics. The AMPT 
model can be extended, e.g., to include hydrodynamic 
evolution at the early stage when local thermalization 
is likely, in order to conveniently study the equation of 
state of the partonic matter. Experiments at RHIC such 
as d-|-A can also help to reduce the uncertainties in the 
physical input such as the parton distribution functions 
in heavy nuclei. With continuing efforts in both theoreti- 
cal and experimental heavy ion physics, we hope that this 
multi-phase transport model will eventually incorporate 
the essential elements of the underlying theory of QCD 
to provide a reliable description of different observables 
in heavy ion collisions within one coherent picture, and 
help us to learn from relativistic heavy ion collisions the 
properties of the quark-gluon plasma formed during the 
early stage of the collisions. 
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APPENDIX: AMPT USERS' GUIDE 

The default AMPT model (current version 1.11) and 
the AMPT model with string melting (current version 
2.11) both use an initialization file 'input. ampt'. The 
analysis directory 'ana/' contains the resulting data files. 
The final particle record file is 'ana/ ampt. dat'. The ver- 
sion number of AMPT is written to both 'ana/version' 
and 'nohup.out' files. The AMPT source code has been 
tested for both f77 and pgf77 compilers on the UNIX, 
Linux, and OSFl operating systems. 

To run the AMPT program, one needs to: 

1. set the initial parameters in 'input. ampt'. If one 
prefers to use run-time random number seed, set 
'ihjsed=ll'. In this way, every run is different even 
with the same 'input. ampt' file. 

2. type 'sh exec &' to compile and run the executable 
'ampt' with some general information written in 
'nohup.out'. 

Key initial parameters in 'input. ampt' are: 

EFRM: y/sWN in GeV, e.g. 200 for the maximum RHIC 
energy. 

NEVNT: the total number of events. 

BMIN, BMAX: the minimum and maximum impact pa- 
rameters (in fm) for all events with BMAX having 
an upper limit of HIPR1(34)-FHIPR1(35) (=19.87 
fm for d-|-Au collisions and 25.60 fm for Au-|-Au 
collisions). 

ISOFT: choice of parton-hadron conversion scenario. 
=1: default AMPT model (version l.a;); 
=4: the AMPT model with string melting (version 
2.2/). 



Note that values of 2, 3, and 5 have never been 
used for publications. They are tests of other string 
melting scenarios: 

=2: a string is decomposed into q+qq-|-minijet par- 
tons instead of using the Lund fragmentation; 
=3: a baryon is decomposed into q-l-qq instead of 
3 quarks; 

=5: same as 4 but partons freeze out according to 
local energy density. 

NTMAX: the number of time-steps for hadron cascade, 
default(D)=150. Note that NTMAX=3 effectively 
turns off hadron cascade, and a larger than the de- 
fault value is usually necessary for observables at 
large rapidity or large pseudorapidity. We use NT- 
MAX=1000 for HBT studies in central Au-hAu col- 
lisions due to the need for the space-time informa- 
tion of last interactions and for LHC calculations 
due to the longer lifetime of the formed matter. 

DT: value of the time-step (in fm/c) for hadron cascade, 
D=0.2. Note that tout = NTMAX x DT is the 
termination time of hadron cascade. 

PARJ(41): parameter a in the Lund symmetric frag- 
mentation function. 

PARJ(42): parameter b in the Lund symmetric frag- 
mentation function (in GeV^^). Note that we use 
defauh value in HIJING (a = 0.5 and b = 0.9) 
for d-|-Au collisions, and a = 2.2 and b = 0.5 for 
collisions of heavy nuclei. 

flag for popcorn mechanism: D=l(Yes) turns on the 
popcorn mechanism. In general, it increases baryon 
stopping. 

PARJ(5): controls BMB vs. BB in the popcorn mecha- 
nism, D=1.0. 

shadowing flag: D=l(Yes) turns on nuclear shadowing. 

quenching flag: D=0(No) turns off jet quenching since 
the parton cascade ZPC simulates final-state ef- 
fects. 

pO cutoff: D=2.0 (in GeV/c) for po in HIJING for mini- 
jet production. 

parton screening mass: controls the parton cross section. 
D=3.2264 (in fm~^). Its square is inversely propor- 
tional to the parton cross section. Use D=3.2264 
for 3mb, and 2.2814 for 6mb. 

ihjsed: choice of the random number seed, D=0. 

=0: take the 'Ran Seed for HIJING' in 'input.ampt' 
and disregard the random value generated in the 
file 'exec'. 

= 11: take the HIJING random seed at runtime 
from the file 'exec', with the seed written in 'no- 
hup.out' and 'ana/version'. 
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Ran Seed for HIJING: random number seed for HIJING 
when ihjsed=0. 

Kshort decay flag: depends on the experimental correc- 
tion procedure, D=0 turns off Kshort decays after 
the hadron cascade. Note that decays of the fol- 
lowing resonances and their antiparticles are always 
included: p, uj, tj, K*, </), A, N*(1440), N*(1535), 
E° (in order to include its feed down to A). 

optional OSCAR output: if set to 1, outputs 
in OSCAR1997A format [13 are written in 
'ana/parton.oscar' and 'ana/hadron.oscar'. 

dpcoal: parton coalescence distance in momentum space 
(in GeV/c). 

drcoal: parton coalescence distance in coordinate space 
(in fm). dpcoal, drcoal both have D=10^ for 
nearest-neighbor coalescence in the AMPT model 
with string melting. 

Key output file are: 

ana/ampt.dat: It contains particle records at hadron 
kinetic freeze-out, i.e., at the last interaction. For 
each event, the first line gives: event number, test 
number(=l), number of particles in the event, im- 
pact parameter, total number of participant nucle- 
ons in projectile, total number of participant nu- 
cleons in target, number of participant nucleons in 
projectile due to elastic collisions, number of par- 
ticipant nucleons in projectile due to inelastic col- 
lisions, and corresponding numbers in target. Note 



that participant nuclcon numbers include nucleons 
participating in both clastic and inelastic collisions. 
Each of the following lines gives: PYTHIA particle 
ID number, three-momentum(pj:,pj,,p2), mass, and 
space-time coordinates(a;,?/,2:,i) of one final particle 
at freeze-out. Note that momenta are in units of 
GeV/c, mass in GeV/c^, space in fm, and time in 
fm/c. If a particle comes from the decay of a reso- 
nance which still exists at the termination time of 
hadron cascade, then its space-time corresponds to 
the decay point of the parent resonance. Also note 
that the x-axis in the AMPT program is defined as 
the direction along the impact parameter, and the 
z-axis is defined as the beam direction. 



ana/zpc.dat: similar to 'ana/ampt.dat' but for partons. 
The first line of each event gives: event number, 
number of partons in the event, impact parame- 
ter, number of participant nucleons in projectile 
due to elastic collisions, number of participant nu- 
cleons in projectile due to inelastic collisions, and 
corresponding numbers in target. Each of the fol- 
lowing lines gives: PYTHIA particle ID number, 
thvee-momentumlpx ,Py,Pz), mass, and space-time 
coordinates(a;,?/,z,t) of one final parton at freeze- 
out. 
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